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Abstract: We use the holographic dual of a finite-temperature, strongly-coupled, gauge 
theory with a small number of flavors of massive fundamental quarks to study meson excita- 
tions and deep inelastic scattering (DIS) in the low-temperature phase, where the mesons are 
stable. We show that a high-energy flavor current with nearly light-like kinematics disappears 
into the plasma by resonantly producing vector mesons in highly excited states. This mech- 
anism generates the same DIS structure functions as in the high temperature phase, where 
mesons are unstable and the current disappears through medium-induced parton branching. 
To establish this picture, we derive analytic results for the meson spectrum, which are exact 
in the case of light-like mesons and which corroborate and complete previous, mostly numer- 
ical, studies in the literature. We find that the meson levels are very finely spaced near the 
light-cone, so that the current can always decay, without a fine-tuning of its kinematics. 
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1. Introduction 



Motivated by some experimental results at RHIC, which suggest that the deconfined matter 
produced in the intermediate stages of a ultrarelativistic nucleus-nucleus collision might be 
strongly interacting, there is currently a large interest towards understanding the properties 
of strongly coupled field theories at finite temperature within the framework of the AdS/CFT 
correspondence (see the review papers |, |, |, g, | for details and more references). A 
substantial part of this effort has been concentrated on studying the response of such a plasma 
to energetic, 'hard', probes, so like heavy quarks M, 
quark-antiquark pairs) 0, 0, |l|, |0|, [2l], [22|, 
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or photons |2g, 27 



mesons (or 

in an attempt to elucidate some intriguing RHIC data, so like the unexpectedly large 'jet 
quenching', or to provide alternative signatures of a strongly-coupled matter. 



In particular, AdS/CFT calculations of deep inelastic scattering (DIS) |32|, |33|, |26|, |27 
1, 1, 1, | 



37, pq] have given access to the structure of strongly coupled matter at high 



- 1 - 



energy and for small space-time separations and thus revealed an interesting picture, which 
is quite different from the corresponding picture in a gauge theory at weak coupling. One has 
thus found that there are no point-like 'partons' at strong coupling, that is, no constituents 
carrying a sizeable fraction x ~ 0(1) of the total longitudinal momentum of a 'hadron' [32. 33] 
or plasma ^7]] at high energy. This has been interpreted as the result of a very efficient 
branching process through which all partons have fallen at the smallest values of x which are 
consistent with energy conservation. This interpretation is consistent with arguments based 
on the operator product expansion at strong coupling [32, 39] and is further supported by 
the fact that the DIS structure functions were found to be large at small values of x <C 1, 
where they admit a natural interpretation in terms of partons [33, [26|, |3|. The central scale 
in this picture is the saturation momentum Q s (x), which defines the borderline between the 
large-x and large virtuality (Q > Q s (x)) domain which is void of partons, and the small-x, 
small virtuality {Q < Q s (x)) domain where parton exist, with occupation numbers of order 1 
(a situation somewhat reminiscent of, but more extreme than, parton saturation in QCD at 
weak coupling pCf| ). This scale Q s plays an essential role also for other high energy processes, 
of direct relevance to heavy ion collisions, like the energy loss and the momentum broadening 
of a heavy quark ]i~5|, [l~6[j. It can furthermore be related to the dissociation length for a large, 



semiclassical, 'meson' [18, 19, EG, 21]. 



So far, most studies of DIS at finite temperature and strong coupling were concerned 
with the M = 4 supersymmetric Yang-Mills (SYM) theory, where the role of the virtual 
photon is played by the 7£-current — a conserved current associated with a U(l) global 
symmetry which couples to massless fields in the adjoint representation of the color group 
SU(iV c ). Very recently, the same problem has been addressed [31] within the context of a 
J\f = 2 supersymmetric plasma obtained by adding iVf hypermultiplets of (generally massive) 
fundamental fields to the M = 4 SYM plasma, in the 'probe' limit where Nf <C N c . In that 
case, the 'photon' is a flavor current which couples to a pair of fields (fermions and scalars) 
in the fundamental representation of SU(iV c ), that we shall refer to as 'quarks'. In order to 
describe the results in Ref. [[H]] and also our subsequent results in this paper, it is useful to 
briefly recall the structure of the holographic dual of the N = 2 plasma at strong coupling 
A = g 2 N c > 1, also known as the 'D3/D7 model' fl], ||, ||, | (see also Sect. | below). 

The supergravity fields (in particular, the Abelian gauge field dual to the flavor current) 
live in the worldvolume of one of the Nf D7-branes that have been inserted in the AdSs x S 5 
Schwarzschild background geometry dual to the M = 4 SYM plasma. For the case where the 
fundamental fields are massive, the D7-branes are separated in the radial direction from the 
iV c D3-branes located at the 'center' of AdSs x S 5 . The distance between the two systems 
of brane then fixes the 'bare' mass of the fundamental 'quarks', represented by Nambu-Goto 
strings stretching from a D7-brane to a D3-brane. Flavorless 'mesons', or quark-antiquark 
bound states, can be described either as strings with both endpoints on a D7-brane, or (at 
least for sufficiently small meson masses and spins) as normal modes of the supergravity fields 
propagating in the worldvolume of the D7-brane. In this paper, we shall adapt the second 
point of view, that of the normal modes. 
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Figure 1: Dispersion relation ui(k) for the lightest spin-zero mesons in the low temperature phase, 
or Minkowski embedding |||/. The solid blue curve corresponds to a pseudo-scalar meson, whereas 
the red dashed curve corresponds to a scalar meson. The solid black line corresponds to lo = k. 



For zero or sufficiently low temperatures, the mesons are strongly bound [43, 23]: their 
binding energy almost compensates the large quark masses (proportional to the string ten- 
sion), so that the meson masses remain finite in the strong coupling limit A — > oo (the 
'supergravity approximation'), to which we shall restrict ourselves. In that limit, the quarks 
become infinitely massive and then the mesons are stable: they form an infinite, discrete, 
tower of (scalar and vector) modes, distinguished by their quantum numbers. A remarkable 
property of the meson spectrum that will play an essential role in what follows is that, at 
finite temperature, the dispersion relation for a given mode changes virtuality, from time-like 
to space-like, with increasing momentum [23|. This is illustrated in Fig. In particular, 



there exists an intermediate value of the momentum at which the mode becomes light -like. 

This situation persists for sufficiently low temperatures, so long as the D7-brane, although 
deformed by the attraction exerted by the black hole, remains separated from the latter. But 
with increasing temperature, one finds a first-order phase transition (from the 'Minkowski 
embedding' to the 'black hole embedding') at some critical temperature T c , at which the tip 



of the D7-brane suddenly jumps into the black hole (BH) horizon [44, 23|. For T > T c , all 
the mesons 'melt' : their dispersion relations acquire large imaginary parts (comparable to 
their real parts) , showing that the bound states are now highly unstable [^4|, |25| . 

One should stress that this peculiar 'meson melting' phase transition is specific to this 
model and has no analog in QCD. The corresponding situation in QCD is not yet fully clear 1 , 



1 For instance, potential models using lattice QCD input predict that all charmonium states and the excited 
bottomonium states dissolve in the QGP juste above T c (see the recent review [0). But the most refined, 
recent, lattice calculations cannot exclude the survival of charmonium states up to T = 1.5T C |46|. 



- 3 - 



but on physical grounds (given that the deconfinement 'phase transition' is truly a cross-over) 
one would expect the mesons to gradually melt when increasing T above T c , according to their 
sizes. The smallest mesons, those built with the heavy quarks, may survive in the temperature 
range pertinent to RHIC or LHC. Assuming the quark-gluon plasma to be effectively strongly 
coupled within this range, one could use the 'low-temperature' phase of this model to get 
some insight into the properties of the still surviving, very small and very heavy, mesons, and 
the 'high-temperature' phase for the larger, lighter, mesons that have already melt. 

Returning to the problem of the DIS for the flavor current, it is quite clear that the 
situation is very different in the high-temperature phase (T > T c ) as compared to low- 
temperature one (T < T c ). At T > T c , the problem is conceptually the same as for the 
7£-current in the N = 4 plasma [26, 27 1. The space-like current fluctuates into a system 
of virtual partons — for the flavor current, this system involves a pair of fundamental fields 
together with arbitrary many N = 4 quanta — whose subsequent evolution depends upon the 
kinematics. If the energy u of the current is low enough (for a fixed virtuality Q), the partonic 
fluctuation closes up again and essentially nothing happens: the space-like current is stable, 
so like in the vacuum, due to energy-momentum conservation. (We ignore here tunnel effects 
at finite temperature, which are exponentially small [26, |36|.) But if the energy is sufficiently 
high, such that to > Q 3 /T 2 , the virtual partons live long enough to feel the interactions with 
the plasma (in the dual gravity problem: the attraction exerted by the black hole) and under 
the influence of these interactions they keep branching until they disappear into the plasma. 
In the dual gravity theory, this medium-induced branching is seen as the fall of the dual gauge 
field into the BH horizon. Since this process is fully driven by interactions with the N = 4 
plasma, the saturation momentum and also the DIS structure functions at small x should be 
the same as for the 7£-current, up to a global factor which counts the number of degrees of 
freedom to which couples the current for the 7£-current and, respectively, N C N{ for the 
flavor one). This is indeed what was found in Ref. [31] for the high-temperature phase. 

But at lower temperatures T < T c , the situation turns out to be more subtle. On the 
supergravity side, it is a priori clear, from geometrical considerations, that the gauge field 
dual to the flavor current cannot fall into the black hole, for any energy: indeed, the support 
of this field is restricted to the worldvolume of the D7-brane, which is now separated from the 
BH horizon along the radial direction. This reasoning led the authors in Ref. [31] to conclude 
that DIS should not be possible in this case, however high is the energy. It is understood 
here that the high-energy limit does not commute with the large-A limit: the energy must 
remain small not only as compared to the quark mass m q ~ \/AM gap , but also as compared 
to the mass of the lowest string excitations m str ing ~ A 1//4 Af gap |43|| , which are not described 
by the supergravity fields. Here, M gap is the lowest meson mass at zero temperature, and is 
independent of A (see Sect. || below). 

Our main, new, observation is that, although it cannot fall into the BH, the flavor current 
can nevertheless disappear into the plasma by resonantly producing space-like vector mesons 
which, as already mentioned, are indeed supported by the plasma. (For a time-like current 
in the vacuum, the resonant production of mesons has been discussed in App. B of Ref. |25[|.) 
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For this process to be possible, the kinematics of the current should match with the dispersion 
relations for the vector mesons. For this interaction to qualify as 'deep inelastic scattering', 
the associated structure functions — determined by the coupling of the current to the mesons 
and computed as the imaginary part of the current -current correlator — must be non-zero 
in a continuous domain of the phase-space, and not only at discrete values of the energy. In 
this paper, we shall demonstrate that these conditions are indeed satisfied at sufficiently high 
energy. Our final result is that the structure functions for flavor DIS are exactly the same 
in this low temperature-phase as in the high-temperature phase, although the respective 
physical pictures are quite different. This result is in fact natural, as we shall later explain. 

To develop our arguments, we shall perform a detailed study of the meson excitations in 
the high-energy, space-like, kinematics relevant for DIS off a strongly coupled plasma; that 
is, uj 3> Q 3> T, with Q 2 = k 2 — uj 2 > 0. We shall focus on vector mesons with transverse 
polarizations, which provide the dominant contribution to DIS at high energy [26|, but we 
expect similar results to apply for other types of excitations (longitudinal vector mesons, 
scalar and pseudoscalar ones). Also, for technical reasons, we shall limit ourselves to the case 
of very heavy mesons, or very low temperature, M gap 3> T, which however captures all the 
salient features of the general situation. 

Concerning the kinematics, we shall find that a space-like current can excite mesons only 
for high enough energies and relatively small virtualities, such that the current and the mesons 
are nearly light-like. This is so because a current with large space-like virtuality encounters 
a potential barrier near the Minkowski boundary (associated with energy-momentum con- 
servation) and thus cannot penetrate in the inner region of the D7-brane, where mesons 
could be created. However, for high enough energy uj > Q 3 /T 2 , this barrier is overcome by 
the gravitational attraction due to the black hole (i.e., by the mechanical work done by the 
plasma 27]), and then the current can penetrate inside the bulk and thus excite mesons. 
The corresponding kinematics being nearly light-like, uj ~ k, we shall focus our attention on 
the respective region of the meson dispersion relation in Fig. |l], but we shall provide analytic 
approximations also for the other regions (in Sect. |4|). Our results are as follows. 

For the strictly light-like mesons (uj = k), we shall construct in Sect. [| exact, analytic, 
solutions for the spectrum and the wavefunctions, which take particularly simple forms for 
large quantum numbers n>l. We shall thus find an infinite tower of equally spaced levels, 
with high energies and a large level spacing: uj n — tiAuj where ~ T(Mg Si p/T^ A^/g a p. 
(For comparison, at zero momentum, the energy of the mode n is uj n (k = 0) ~ nM gap .) 
Similarly, for the gauge field dual to a light-like flavor current, we shall find exact 'non- 
normalizable' solutions, from which we shall compute the retarded current-current correlator 
in the high energy limit uj 3> T(M gap /T) 3 . As expected, this propagator exhibits poles at 
the energies of the light -like mesons, so its imaginary part is an infinite sum over delta- 
like resonances. The coefficient of each delta-function represents the probability for the 
resonant production of a meson by a current whose energy is exactly u n . Conversely, they 
also describe the rate for the decay of a vector meson into an on-shell photon, a mechanism 
recently proposed as a possible signal of strong coupling behaviour in heavy ion collisions [3C]. 
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The resonant production of mesons remains possible also for slightly space-like kinemat- 
ics, and, of course, for any time-like kinematics, but this is perhaps not the most interesting 
physical situation, as it requires the energy of the current to be finely tuned to that of a meson 
mode. Given the large level spacing Aw indicated above, it looks at a first sight unlikely that 
a small uncertainty 5uj <C u in the energy of the current — as inherent in any scattering 
experiment (even a Gedanken one !), where the 'photon' is not a plane wave, but a wave 
packet — could help reducing the need for the fine-tuning. If that was true, it would mean 
that for the whole high-energy region in phase-space, except for a set of zero measure (as de- 
fined by the dispersion relations for the meson modes), the current survives in the plasma for 
arbitrarily long time. However, as we shall argue now, that conclusion would be a bit naive, 
as it underestimates the consequences of a small fluctuation in the energy, or the virtuality, 
of the current for the problem at hand. 

The main point is that the energy uncertainty 5uj should not be compared to the (rela- 
tively large) level spacing Au between two successive resonances, but rather to the change 
in energy which is necessary to cross from one meson level to another at a fixed value k of 
the momentum. Indeed, one should not forget that the dispersion relation in the relevant 
kinematics is nearly light-like, that is, uj n ~ k n for the nth mode. Hence, when increasing the 
energy by Au; to move from one level to a neighboring one, one is simultaneously increasing 
the momentum k by the same, large, amount — one moves along the light-cone. But in a 
scattering problem, the momentum k of the current is fixed and its energy u> has generally 
an uncertainty 8u> , related to the fact that the source producing the current has been acting 
over a finite period of time 5t: 5uj ~ 1/St. Before we discuss this time 5t, let us make the 
crucial observation that the meson levels are very finely spaced in energy when probed at a 
fixed value of k. This is a general feature of the high-energy kinematics, which is further 
amplified by the peculiar shape of the meson dispersion relation near the light -cone. 

As a simpler example, recall first the situation at zero temperature [43], where the meson 
dispersion relation reads, schematically, u) n (k) = (k 2 + n 2 Mg ap ) 1 / 2 ~ k + n 2 Mg ap /2/c, with 
the approximate equality holding when k S> nAf gap . Hence the energy jump 5uj n {k) = 
uj n+ \(k) — uj n (k) needed to cross from one mode to another at fixed k is 5uj n {k) ~ nM 2 ap /fc 
and becomes smaller and smaller when increasing k. As anticipated, the modes are very finely 
spaced at large k. Returning to the finite-T case of interest, it turns out that the respective 
dispersion relation is even more sensitive to small variations in the virtuality of the mode near 
the light-cone. Specifically, we shall find in Sects. [D^ and |6| that the level spacing at fixed k 
defined as above scales with k as 5co n (k) ~ T{T/k) 1 ^ when k ~ k n = nAu (see Fig. |J). As 
anticipated, this is considerably smaller than the energy spacing Au ~ r(M gap /r) 3 at fixed 
virtuality : indeed, Au)/8u) n ~ n 1 / 3 (M gap /T) 4 » 1. 

To understand the typical energy uncertainty 5lo of the current, one needs an estimate 
for its interaction time in the plasma t- in t. Indeed, the source producing the current should 
act over a comparatively shorter time 5t < t- m t in order for the subsequent dynamics to be 
observable. Via a time-dependent analysis of the dynamics of the dual gauge field in Sect. ^, 
we shall find that ii n t is controlled by the progression of the gauge field within the D7-brane, 
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Figure 2: Qualitative illustration of the meson dispersion relation uj(k) in the vicinity of the light- 
cone. Four successive modes are shown. Although widely spaced along the light-cone, the modes are 
close to each other at any given value of k. Thus, in order to jump from one mode to a neighboring 
one, one needs a relatively large energy jump Alo at fixed virtuality Q, but only a small energy jump 
5w at fixed momentum k. 



which yields t; n t ~ (/c/T) 1//3 /T (similarly to the 7£-current [^]]). This estimate implies a lower 
limit 5co > 1 /ti n t on the energy uncertainty of the current which is of the order of the level 
spacing 5uj n indicated above. This justifies performing an average over neighboring levels 
in the calculation of the imaginary part of the current-current correlator. This averaging 
smears out the meson resonances and produces our main result in this paper, Eq. ( 5.1§| ). As 
anticipated, this result is identical to the DIS structure functions in the high-temperature 
phase, which shows that the current is completely absorbed by the plasma in both cases. 

The analysis in Sect, ^jalso allows us to deduce a space-time picture for the nearly light- 
like mesons in the semiclassical regime at large quantum numbers n>l, where the notion of a 
classical orbit makes sense. We thus find that the period for one orbit is At n ~ (tj n /T) 1 / 3 /T, 
where we recall that the energy of the bound state is u n = nAco. Furthermore, we find that 
the meson spends the major part of this time far away from the tip of the D7-brane, at 

1/3 

relatively large radial distances ~ u> n . This is so because its orbital velocity is much higher 
near the tip than at larger radial distances. It is finally interesting to notice that, in this 
light-like kinematics, the period At n of the bound state has the same parametric dependence 
upon its energy as the interaction time t- m t of the current, and similarly for the typical radial 
location of the meson versus the saturation momentum Q s (k) ~ k 1 / 3 for the current. 
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2. Mesons in the D3/D7 brane model at finite temperature 



According to the AdS/CFT correspondance [^7], 48, 49 1, the four-dimensional J\f = 4 super- 
Yang-Mills (SYM) gauge theory with 'color' gauge group SU(iV c ) is dual to a type IIB 
string theory living in the ten-dimensional curved space-time AdSs x S 5 , which describes the 
decoupling limit of N c black D3-branes. By further adding a black brane to this geometry, 
one obtains the holographic dual of the finite-temperature, plasma, phase of the M = 4 SYM 
[pCf| . The ensuing metric reads (see, e.g., Q) 

ds 2 = (-f(u)dt 2 + dx 2 ) + ^(^- + u 2 dnl) , (2.1) 



L 2 v J v ' ' > u 2 \ f(u) 

where f(u) = 1 — Uq/u 4 , with no = ttL 2 T the radial position of the black hole horizon 
and T the common temperature of the N = 4 SYM plasma and of the black hole. The 
curvature radius L is defined in terms of the string coupling constant g s and the string length 
scale l s via L 4 = 47rg s N c The holographic dictionary relates the gauge and string theory 
coupling constants as g 2 = Airg s . In the "strong coupling limit" of the gauge theory, defined 
as N c — > oo, A = g 2 N c — ► oo, with g fixed and small (g <C 1), the string theory reduces to 
classical supergravity theory in the AdSs x S 5 Schwarzschild geometry with metric fl2.1|). 

All fields in the N = 4 SYM theory transform in the adjoint representation of SU(iV c ). 
Fields transforming in the fundamental representation of the gauge group can be introduced in 
the gravity dual by inserting a second set of D-branes in the supergravity background [^, f|. 
In particular, we consider the decoupling limit of the intersection of N c black D3-branes and 
N{ D7-branes as described by the array: 

0123456789 
D3 x x x x (2.2) 

D7 xxxxxxxx 

where the first four dimensions (0, 1, 2, and 3) correspond to the Minkowski coordinates 
{t,x 1 } and the last six ones (from 5 to 9) to the six-dimensional space with coordinates 
{u, ris}. The dual field theory is now an N = 2 gauge theory consisting of the original SYM 
theory coupled to Nf fundamental hypermultiplets which consists of two Weyl fermions and 



their superpartner, complex, scalars (see e.g. [51]). For brevity, we shall globally refer to 
these fundamental fields as 'quarks'. In the limit where the number of flavors is relatively 
small, N{ -C N c , the D7-branes may be treated as probes in the black D3-brane geometry 
d2.1|) . That is, the D7-branes are generally deformed by their gravitational interactions with 
the D3-branes and the black hole, but one can neglect their back reaction on the ambient 



geometry, Eq. (2.1). The ensuing geometry is dual to a N = 2 plasma at finite temperature in 
which the effects of the fundamental degrees of freedom (say, on thermodynamical quantities) 
represent only small corrections, of relative order g 2 N{ = \(N{/N C ) <C 1 (see e.g. |23|] ). 

Although both the D3-branes and the D7 ones fill the Minkowski space, these two types of 
branes need not overlap with each other, as they can be separated in the 89-directions, which 
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Figure 3: Schematic representation of the geometry of the D3-D7 model at finite temperature. This 
drawing is adapted to the low-temperature, or very heavy quark, situation, where Rq 3> po — 1 ■ See 
the text for more explanations. 

are orthogonal to both of them. When this happens, the conformal symmetry is explicitly 
broken already at classical level 2 and then the fundamental fields in the dual gauge theory 
become massive: their 'bare' mass is proportional to the radial separation u m between the 
two sets of branes at zero temperature. Indeed, a fundamental field is 'dual' to an open string 
connecting a D7-brane to a D3-brane, so its 'bare' mass is equal to the string length u m 
times the string tension: 

mq= 2^ S =VX 2^- (2 ' 3) 
To render such geometrical considerations more suggestive, it is helpful to perform some 



changes of coordinates [44, 23]. First, we introduce a new, dimensionless, radial coordinate 
p, related to the coordinate u via 3 



{uqpY = u z + Ju 4 - u A Q . (2.4) 



2 Quantum mechanically, the conformal symmetry is broken by the D7-branes even when they overlap with 
the D3-branes, i.e., when u m = 0. But the /3-function for the 't Hooft coupling A = g 2 N c is of order Nf/N c 
and thus is suppressed in the probe limit Ni/N c — > 0. 

3 We notice that p is related to the Fefferman-Graham ]k| radial coordinate z via zjsfl = L 2 /(uop). 
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Note that that the BH horizon corresponds to po = 1 and the Minkowski boundary to p — > oo, 
with uqp ~ \plu when u^> uq (i.e., p S> 1). Then the background metric (2.1) becomes 



f 2 

■^dt 2 + fdx 2 
f 



^[dp 2 + p 2 dO|] , (2.5) 



where 



f(p) = l-y 4 , f(p) = l + y i . (2.6) 

It is furthermore useful to adapt the metric on the five-sphere to the D7-brane embedding. 
Since the D7-brane spans the 4567-directions, we introduce spherical coordinates {r, f^} in 
this space and {R, 4>} in the orthogonal 89-directions. Denoting by 9 the angle between these 
two spaces, we have (see also Fig. 0) 



and therefore 



p 2 = r 2 + R 2 , r = pcos6, R = psin6, (2.7) 

dp 2 + p 2 dnj = dp 2 + p 2 (d6 2 + cos 2 9 dQ% + sin 2 9 d0 2 ) 

= dr 2 + r 2 d^l + di? 2 + R 2 d<p 2 . (2.8) 

Note that, on the D7-brane, the Minkowski boundary lies at r — ► oo. 

To specify the D7-brane (background) embedding, we require translational symmetry in 
the 0123-space and rotational symmetry in the 4567-directions, and fix <fi = 0. Then the 
embedding can be described as the profile function R = R v (r). The subscript V on R v 
stays for the 'meson vacuum': the small fluctuations of the D7-brane around its stationary 
geometry are dual to low-lying 'mesons' in the boundary gauge theory, i.e., (colorless and 
flavorless) bound states which involve a pair of fields from a fundamental hypermultiplet 
- say, a quark-antiquark pair. Such mesons are represented by strings with both ends on 
the D7-branes and thus can be studied (at least for small enough meson sizes and masses; 
see below) by examining the small fluctuations of the worldvolume fields on the D7-branes. 
These include the fluctuations 5(f) and 5R in the shape of the D7-brane — which give rise to 
pseudo-scalar and scalar mesons, respectively — , and also fluctuations of the worldvolume 
gauge fields, which describe vector mesons. The 'vacuum' profile function R v (r) and the 
spectrum of the various type of fluctuations have been systematically studied in the literature, 
via analytic methods in the zero-temperature case jE]], and via mostly numerical methods 
at non-zero temperature [42, 44, 23, p3, pH]. In what follows, we shall collect the previous 



results which are relevant for the present analysis, with a minimum of formulas. 

The dynamics of the D7-brane is described the Dirac-Born-Infeld (DBI) action 4 . The 
profile function R v (r) for the 'vacuum' embedding is obtained by solving the equation of 



4 The full D7-brane action also involves a Wess-Zumino term, but this plays a role only for those gauge 
field configurations having non-trivial components along the three-sphere S 3 internal to the D7-brane |f3| , p5| . 
Such fields do not enter the study of deep inelastic scattering and will not be considered throughout this work. 
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motion for R(r) which follows from this action. The meson spectrum is then obtained by 
solving the linearized equations of motion (EOM) which follow after expanding the DBI action 
to quadratic order in small fluctuations around the 'vacuum' embedding. 

At zero temperature, one finds that the 'vacuum' profile is trivial, i.e., independent of r : 

Rv{f) = Rq = V2 — — (zero temperature) , (2-9) 

where we recall that u m is the separation between the two types of brane in the original 
radial coordinate u. (At T = 0, Eq. ( ]2.4| ) reduces to p = \/2{u/uq) where uq is an arbitrary 
reference scale, which drops out from the final results.) The EOM for the small fluctuations 
have been solved exactly, in terms of hyper geometric functions [|43|]. At zero temperature, 
both Lorentz symmetry and supersymmetry are manifest. Accordingly, for a meson with four- 
momentum = (u, 0, 0, k), the dispersion relation u(k) involves only the 'invariant mass' 
combination M 2 = lo 2 — k 2 . Besides, this relation depends upon two 'quantum numbers': 
a 'radial' number n = 0,1,2,..., which counts the number of zeroes of the corresponding 
wavefunction in the interval < r < 00, and an 'angular' number £, with t = 0,1,2,..., 
which refers to rotations along the S 3 component of the D7-brane. (In the dual gauge 
theory, I represents a charge under the internal symmetry group SO (4) which is dual to 
rotations on S 3 .) Supersymmetry together with the global SO(4) symmetry imply additional 
degeneracies for the meson spectrum, as discussed in Specifically, Ref. |43| found the 

following dispersion relation 

M 2 (n,£) = u 2 ni {k)-k 2 = ^A( n + £ + l)(n + £ + 2) (2.10) 

for both (pseudo)scalar and vector mesons. Note the presence of a mass gap in the spectrum: 
the mass of the lightest mesons is non-zero, namely, 

M gap = 2^2^. (2.11) 

Note also that the meson masses are much smaller, by a factor l/VX, than the bare quark 
mass, Eq. (|2.3|). This shows that in this strong coupling limit the mesons are tightly bound: 
in the total energy, the binding energy almost cancels the mass of the quarks. 

At finite temperature, the D7-brane feels the attraction exerted by the black hole and thus 
is deflected towards the latter — the stronger the deviation, the shorter is the radial separation 
u (or p) between the two. This deflection becomes negligible towards the Minkowski boundary 
(r — > 00), where the profile function R v (r) approaches the value Rq that it would have (at 



any r) at T = 0. More precisely, for asymptotically large r one finds E3, 44, u. 



R v (r) ~ R - 4 , (2-12) 



with Ro related to the 'bare' quark mass, as in Eq. (|2.3|) , and c a positive number proportional 
to the quark condensate. 
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On the other hand, closer to the black hole horizon (p ~ 1), one finds two different types 
of behaviour — corresponding to two thermodynamically distinct phases separated by a first- 
order phase transition — , depending upon the ratio M gap /T = 2ttRq between the (zero-T) 
mass gap and the temperature: (i) for relatively large values of i?o, larger than a critical 
value numerically found as R c — 1.306 [ffijfl , the D7-branes close off above the black hole 
horizon ('low-temperature', or 'Minkowski embeddings'); (ii) for Rq < R c , the D7-branes 
extend through the horizon ('high-temperature', or 'black hole embeddings') 5 . In the gauge 



theory, the most striking feature of this transition is the change in the meson spectrum [44] : 



in the low temperature phase, the spectrum of mesons has a mass gap and the bound states 



are stable, so like at T = [23|; in the high temperature phase, there is no mass gap and 
the mesonic excitations are unstable and characterized by a discrete spectrum of quasinormal 
modes (i.e., they have dispersion relations with non-zero, and large, imaginary parts) [f24|, |25f] . 

For the reasons explained in the Introduction, in this paper we shall restrict ourselves to 
the low-temperature phase, in which the mesons are stable. The corresponding dispersion 



relations have been numerically computed in Ref. [23|, at least within restricted regions of 
the phase space. As expected, the spectrum shows deviations from both Lorentz symmetry 
and supersymmetry, and these deviations become more and more important with increasing 
temperature (for a given M gap ). What was perhaps less expected and, in any case, remarkable 
is the pattern of the violation of the Lorentz symmetry by the spectrum: when increasing the 
momentum A; of a given mode (i.e., , for fixed values of n and £, which remain good 'quantum 
numbers' also at finite temperature), the 'virtuality' — Q 2 = uj 2 Jk) — k 2 of that mode is 
continuously decreasing, from time-like values (-Q 2 > 0) at relatively low k to space-like 
values (— Q 2 < 0) for sufficiently large k, in such a way that, for asymptotically large k, the 
dispersion relation approaches a limiting velocity which is strictly smaller than one: 

u(k) ~ vq k with vq < 1 as k — > oo . (2-13) 



It has been furthermore noticed in the numerical analysis in Ref. |23j that, with increasing 
k, the mode wavefunction becomes more and more peaked near the bottom (r — ► 0) of the 
D7-brane. This led to the interesting suggestion, which was furthermore confirmed by the 
numerical results, that the limiting velocity vq coincides with the local velocity of light at 
r « 0, i.e., p « Rv(0) ■ 



9tt 

9zz 



r=0 



(2.i4) 

f(R,(o)) 



As we shall later argue in Sect. |3|, this identification follows indeed from the respective EOM. 
Our analytic study will also clarify other aspects of the dispersion relation, like the precise 



conditions for the onset of the asymptotic behaviour (2.13) and the subleading corrections to 



5 There is a discontinuous jump in between these two phases: with increasing temperature, the 'tip' of the 
D7-brane which in the Minkowski embeddings lies at R v (0) reaches an absolute minimum value R v (0) ~ 1.15 
corresponding to Ro = R c and then jumps into the black hole horizon; see [E3| for details. 
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it, which in particular contain the dependence upon the quantum numbers n and t. More 
generally, we shall be able to construct piecewise analytic approximations for the dispersion 
relation oo n t(k) and also for the wavefunctions of the modes, which will confirm the numerical 
findings in Ref. [^3| and provide further, analytic, insight into these results. Although, in our 
analysis, we shall cover all kinematical domains in k and thus provide a global picture for the 
meson spectrum, our main focus will be on the nearly light-light mesons with uj ~ k. Indeed, 
as we shall explain in the next section, this regime is the only one to be relevant for the deep 
inelastic scattering of the flavor current. 

3. Deep inelastic scattering off the J\f = 2 plasma 

The N = 2 theory with Nf flavors of equal mass has a global U(Nf) ~ SU(N{) x £7(l) q sym- 
metry (describing flavor rotations of the fields in the fundamental hypermultiplets), to which 
one can associate Nj conserved currents bilinear in the 'quark' operators (see Appendix A in 
|^| for explicit expressions). In particular, the current Jq corresponding to the diagonal sub- 
group U(l)q is associated with the conservation of the net 'quark' number (i.e., the number 
of fundamental quarks and scalars minus the number of antiquarks and hermitean conjugate 
scalars). By adding to the theory a U(l) e . m . gauge field minimally coupled to this Jq 
current (with an 'electromagnetic' coupling which is arbitrarily small), one can construct a 
model for the electromagnetic interactions and thus set up a Gedanken deep inelastic scatter- 
ing experiment which measures the distribution of the fundamental fields inside the plasma. 
One can visualise this process as the exchange of a virtual, space-like, 'photon' (as described 
by the field ^4^) between the strongly coupled N = 2 plasma at finite temperature and a hard 
lepton propagating through the plasma. 

3.1 Equations of motion in the D3/D7 brane model 

Within the D3/D7 brane model, the flavor current Jq is dual to an abelian gauge field A m 
living in the worldvolume of the D7-brane, whose dynamics is encoded in the DBI action. 
According to the gauge/gravity duality, the correlation functions of the operator Jq are 
obtained from the 'non-renormalizable' modes of the field A m , that is, the solutions to the 
classical EOM in the bulk of the D7-brane which obey non-trivial (Dirichlet) boundary 
conditions at the Minkowski boundary: as r — * oo, the solution A m must approach the 
U(l) e .m. gauge field A^ which acts as a source for the current Jq. This should be contrasted 
to the 'normalizable' modes dual to vector mesons, which must vanish sufficiently fast when 
approaching the Minkowski boundary (see below for details). 

In particular, the DIS cross-sections (or 'structure functions') are obtained from the 
(retarded) current-current correlator 

IT-(g) ee i f dV-^^oKfJqXy), J£(0)])t, 

where the brackets (• • • )t denote the thermal expectation value in the M = 2 plasma. To 
compute this two-point function, it is enough to study the linearized EOM for the bulk field 
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A m , i.e., the same equations which determine the spectrum of the low-lying vector mesons, 
but with different boundary conditions at r — > oo. 

Specifically, the polarization tensor (3.1) can be given the following tensorial decomposi- 
tion (in a generic frame) : 

IV(g,T)= L^-^\n x { X ,Q 2 ) + L^-qJ^\ (^-^^II 2 (x,Q 2 ), (3.2) 

where IIi and H2 are scalar functions, is the four-velocity of the plasma in the considered 
frame, Q 2 = q^q^ > is the (space-like) virtuality of the current, and 

Q2 (3.3) 



-2(q-n)T 



is the Bjorken variable for DIS off the plasma. Via the optical theorem, the DIS structure 
functions are obtained as 

Fi(x,Q 2 ) = ^Imll!, F 2 (x,Q 2 ) = ^_^i m n 2 . (3.4) 

In what follows, it will be convenient to compute the (boost-invariant) structure functions by 
working in the plasma rest frame, where = (1,0,0,0) and q^ = (u,0, 0, k), and therefore 
Q 2 = k 2 — to 2 and x = Q 2 /2luT. However, one should keep in mind that the physical 
interpretation of the results is most transparent in the plasma 'infinite momentum frame', 
i.e., a frame in which the plasma is boosted at a large Lorentz factor 7 > 1. Then, the 
kinematic invariants Q 2 and x specify the transverse area (~ 1/Q 2 ) and, respectively, the 
longitudinal momentum fraction (equal to x) of the plasma constituent ('parton') which has 
absorbed the space-like 'photon', and the structure functions represent parton distributions. 
The piece of the DBI action which is quadratic in the gauge fields reads (see e.g. p5|]) 

S 8 = -Pj^T-DrNi J d 8 a^g m Pg^F mn F pq , (3.5) 

where Td7 = 2-K/(2ir£ s ) 8 g s , the space-time indices m,n,p,q run over the eight directions in 
the worldvolume of the D7-brane, g mn is the induced metric on the D7-brane, and F mn = 
d m A n — d n A m . As already mentioned, the EOM must be solved with the following boundary 
conditions 

A m (t, x, r - 00) - A^t, x) = 4°) e -^ +ifcz , (3.6) 

which together with the fact that the equations are linear and homogeneous in all the world- 
volume directions but r imply that the solution A m is such that A r = A S 3 = (i.e., the 
radial and 5' 3 -components of the gauge field are identically zero) and the remaining, four, 
components A^(t, x, r), with [i = t,x,y,z, are plane-wave in the Minkowski directions with 
r-dependent coefficients. Since the gauge fields are independent of the coordinates on S 3 , 
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one can reduce Eq. (|3.5| ) to an effective action in the relevant five dimensions. The induced 
metric in these directions, that we denote as g mn , follows from Eqs. (|2.5|)-(^8|) as 



2,.,., _ 1 ^uppy 



ds 2 (g) 



f 2 

4- dt 2 + fdx 2 
f 



+ ^ (l + Rl) dr 2 , 



(3.7) 



where p 2 = r 2 + R 2 {r), R v = dR v /dr, and we recall that R v (r) is the profile of the D7-brane 
embedding. After integrating over the coordinates on S* 3 , the action (|3.5D reduces to 



S 



{2tt£ 



2\2 



-n 3 T D7 N { / dtd 3 xdr 



Cff( r ) 



9 ~mp~nqp p 
iJ y 1 mn 1 pq > 



(3.8) 



where ^3 = 2tt 2 , m,n,- ■ ■ = t,x,y,z,r, and 



( —X 



Jl . 



9 2 eS {r) \Pv) 



Clearly, the EOM generated by the action ( |3.7| ) read 



gmp-nqp 



PI 



0. 



(3.9) 



(3.10) 



The propagation of the virtual photon along the z axis introduces an anisotropy axis in 
the problem, so the equations of motion look different for the longitudinal (/z = t, z) and 
respectively transverse (p = x, y) components of the gauge field. In what follows we shall 
focus on the transverse fields, Ai with i = x,y, since from the experience with the "/^-current 



27] we expect these fields to provide the dominant contributions to the structure functions 



F\ 2 in the high energy limit. Moreover, our final argument in Sect. || will allow us to also 
reconstruct the flavor longitudinal structure function Fl = F2 — 2xF\ from the corresponding 
one for the 7£-current. The relevant components of the field strength tensor are F r i = d r Ai, 
F t i = dtAi — > —iujAi, and F z i = d z Ai — > ikAi, and the equation satisfied by Aj(r) reads 



Ai + 



d r In 



5cfr 



-9~ g rr~i 



9 



9' 



f 2 



k 



0. 



(3.11) 



This equation must be solved with the Dirichelet boundary condition (3.6) at r — > 00 together 
with the condition of regularity at r = 0. The solution to this boundary- value problem is 
a 'non-normalizable' mode, as opposed to the 'normalizable' modes which describe vector 



meson excitations of the D7-brane: the latter are the solutions Eq. (3.11) which vanish 
sufficiently fast (namely, like A^{r) ~ 1/r 2 ) when r — > 00 |23| . 

Once the 'non-normalizable' solution is known as a (linear) function of the boundary 



value A^ , the current-current correlator (3.1) is obtained, roughly speaking, by taking the 



second derivative of the classical action {i.e., the action ( |3.8| ) evaluated with that particular 
solution) with respect to A\ . This procedure is unambiguous in so far as the euclidean {i.e., 
imaginary-time) correlators are concerned, but it misses the imaginary part for the real-time 
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correlators. Rather, the correct prescription for computing the retarded polarization tensor 
(ETlD reads (for the case of U xx = U yy = U 1 ) g§ 



Ui{x,Q z ) = 



3 d r Aj(r,u;,k) 
Ai(r,uj,k) 



(3.12) 



where i is either x or y. The overall normalization factor reflects the fact that the flavor 
current couples to N C N{ fundamental fields. When using the above formula, the precise 
normalization at r — ► oo, i.e., the boundary value A^' , becomes irrelevant, as it cancels in 
the ratio. Note that, in order to make use of Eq. ( |3.12| ), it is enough to know the solution 
in the vicinity of the the Minkowski boundary. But to that aim, one generally needs to solve 
the EOM for arbitrary values of r, since the second boundary condition is imposed at r = 0. 

In general, the coefficients in Eq. ( [3,llD are rather complicated functions, as visible on 
Eqs. ( |3.7| ) and ( |3.9| ), and this complication hinders the search for analytic solutions. However, 
how we now explain, they can be considerably simplified without loosing any salient feature 
by restricting ourselves to the very low temperature, or very heavy meson, case Rq S> 1, or 
-^gap 3* T. This restriction entails two important types of simplifications. The first one refers 
to the 'vacuum' profile R v (r), which in the general case is known only numerically ^3|, but 
which becomes essentially flat when J?o > 1. Indeed, in that case, the maximal deviation 
from the asymptotic value Rq, namely (see App. A in Ref. p3fl), 

R v (0)-Ro ^ ~2^7 < ( 3 - 13 ) 

is truly negligible, so one can use R v {r) ~ Rq (and hence R v = 0) at any r. 

The second type of simplifications refer to the BH horizon at po = 1 : when Rq ^> 1, 
the condition p 3> 1 is automatically satisfied at any point within the worldvolume of the 
D7-brane. Then, the thermal effects encoded in / and /, which scale like 1/p 4 , cf. Eq. (|2.6| ), 
can be safely neglected in all the terms in Eq. ( |3.11|) except for the last one: indeed, within 
that term, the finite-T deviations 1 — / and 1 — / are potentially amplified by the large energy 
factor oj. Specifically (with p 2 = Rq + r 2 ) 



/ 2 i2 i 2 i 2\ 4/ p 2 r\2 , ^ 2 



where we have also used Q 2 = k 2 — oj 2 . 

To summarize, under the assumption that Rq S> 1, the EOM for the transverse gauge 
fields Ai{r) takes a particularly simple form: 

l; + ~i; + (-^ + ^U = > ( 3 - 15 ) 
r \ P P J 

where Ai = d^4j/dr, p 2 = Rq + r 2 , and we have introduced the dimensionless variables 

oj = uj— = -^=, k = k— = -^, Q 2 = k 2 -oj 2 . (3.16) 



- 16 - 



One should emphasize here that this condition Rq 3> 1 introduces no loss of generality, neither 
for a study of the DIS process (in which case we are anyway interested in u, Q 3> T, and then 
the dominant dynamics takes place at large radial distances p ^ 1 [pq, ^7]]), nor for that of 
the meson spectrum (for which we shall find results which are consistent with the numerical 
analysis in Ref. ||23|| , although that analysis was performed for Rq ~ 0(1)). 

Although considerably simpler than the original equation ( 3.1 1) ) , the above equation is 



still too complicated to be solved exactly, except in the special case Q = 0, to be discussed in 
Sect. |5[ For more general situations, related to either the meson spectrum or the problem of 
DIS, we shall later construct analytic approximations. In preparation for that and in order to 
gain more insight into the role of the various terms in Eq. ( p.15 ), it is useful to first consider 



a different but related problem, whose solution is already known : this is the DIS of the 



^-current 3^, B3L Eg, 27]. 



3.2 Some lessons from the 1Z current 

The "^.-current is a conserved current associated with one of the U(l) subgroups of a global 
SU(4) symmetry of the M = 4 SYM theory. The respective operator is bilinear in the mass- 
less, adjoint, fields of M = 4, and remains conserved even in the presence of the fundamental 
hypermultiplets (i.e., in M = 2 theory), because of the probe limit g 2 Nf <C 1. The super- 
gravity field dual to the 7£-current is, once again, a gauge field A^, whose dynamics however 
is not anymore restricted to the worldvolume of the D7-brane — rather, this field can propa- 
gate everywhere in the AdSs x S" 5 Schwarzschild space-time, in particular, it can fall into the 
black hole. Because of that, the D7-brane plays no role in the case of the 7£-current, so the 
following discussion applies to both J\f = 4 and M = 2 theories (with iVf <C N c , of course). 

For a space-like 7?.-current with high virtuality Q 3> T and for large radial coordinates 
p ^ Pq, the dynamics of the dual 7£-field Ai(r) is described by an equation similar to 



Eq. ( 3.15 ), but where the variables p and r are now identified with each other (since Rq plays 



no role in this case). That is, 

Ai + - Ai + f-^- + Ai = (^-current) , (3.17) 

p V P P H ) 

where now Ai = dAi/dp and it is understood that p> 1. 

The dynamics is driven by the competition between the two terms inside the brackets in 
Eq. ( |3.17| ). The first term, proportional to Q 2 , acts as a potential barrier which opposes to 
the progression of the field towards the interior of AdSs : by itself, this would confine the 
field near the Minkowski boundary, at large radial distances p > pg = Q. The second term, 
proportional to co 2 , is present only at finite temperature (as manifest from its derivation in 



Eq. ( 3.14 )) and it represents the gravitational attraction between the gauge field and the BH. 
For sufficiently small values of p, smaller than p c = (2u)/Q) 1 / 2 , this attraction overcomes the 
repulsive barrier oc Q 2 , and then the overall potential becomes attractive. However, unless 
the energy u is high enough, this change in the potential has no dynamical consequences 6 , 



6 We here ignore the tunneling phenomenon, which is exponentially suppressed |2(], ^] . 
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because the field is anyway stuck near the Minkowski boundary and thus cannot feel the 
attraction. Clearly a change in the dynamics will occur when the energy is so high that 
Pc Pq, which requires u> > Q 3 , or, in physical units, ui > Q 3 /T 2 . When this happens, the 
potential barrier at p — > oo cannot prevent the gauge field to penetrate (through diffusion; 
see the discussion in Sect. |6|) down to the attractive part of the potential at p < p c , and from 
that point on, the potential barrier plays no role anymore. Hence, the dynamics at radial 
distances 1 <C p < p c is controlled by the even simpler equation 

3 ■ 8O 2 

Ai + -Ai H 5- Ai = (fa-current & high energy) , (3.18) 

9 P 

which does not involve Q and thus is formally the same as the equation describing a light- 
light current. This equation can be exactly solved in terms of Airy functions. The general 
solution reads (up to an overall normalization, which is irrelevant) 

Mp) = cAi(-0 + Bi(-0 , e = , (3-19) 

valid for (Q/u) 1 / 3 ) < £ <C cj 2 / 3 . This involves one unknown coefficient c which can be fixed, 
in principle, by matching onto the corresponding solution at smaller distances p ~ 0(1), 
which in particular obeys the appropriate boundary condition at p = po = 1. This boundary 
condition is rather clear on physical grounds: the gauge field can be only absorbed by the 



BH, but not also reflected, hence the near-horizon solution must be a in] "ailing wave |53|, g], 
i.e., a field which with increasing time approaches the horizon. 

The solution near p = 1 obeying this boundary condition can be explicitly computed, 
and its matching onto Eq. ( 3.1SQ can indeed be done [26], but it turns out that this actually 



not needed for the purpose of computing the DIS structure function: the problem of the fa- 
current offers an important simplification, which is worth emphasizing here, since the same 
simplification appears for the flavor current in the high-temperature case (the 'black hole 
embedding') |jT], but not also in the low-temperature, or 'Minkowski', embedding of interest 
for us here. Namely, the infalling boundary condition can be enforced not only near the BH 
horizon, but also at much larger values of p, where Eq. ( 3.1 9| ) applies. This is so since there 



is no qualitative change in the shape of the potential at any intermediate point in the range 
1 < p <C p c which could give rise to a reflected wave. 



The last observation allows us to identify c = i in Eq. ( 3.19 ): indeed, consider this 
approximate solution for p <C u) 1 / 3 , or £ 3> 1, where one can resort on the asymptotic 
expansions for the Airy functions. Using Eq. ( |3.19| ) with c = i, one obtains 

Ai(t,z,p) ~ ^1/4 ex P | - iujt + ikz + i^ £ 3/2 + ijj for 1 < p < cj 1 / 3 , (3.20) 

which is indeed an infalling wave. 

Now that the coefficient c has been fixed, one can use the approximate solution (|3.19[) for 
relatively small values of £ and compute the current-current correlator according to Eq. ( [3.12 ) . 
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(One can adapt Eq. ( |3.12| ) to the 7£-current by multiplying its r.h.s by a factor N C /4:N{.) 
Specifically, Eq. ( p. 19 ) is still correct for £ ~ Q/ui 1 / 3 <C 1, where one can use the small-£ 



expansions for the Airy functions, and thus deduce 

F ^^k&' 3 ' (3.21) 

where in the last estimate we indicated the parametric dependencies of the structure functions 
upon the variables relevant for DIS 7 . As it should be clear from the previous discussion, these 
results hold for sufficiently high energy, u > Q 5 /T 2 , a condition which can be rewritten in 
terms of the DIS variables x and Q 2 as 

- Q ~ < x s {Q) = J, or Q < Q s (x) = -. (3.22) 



2loT ~ Q ' ^ ~ ' x 

On the other hand, for larger values of Bjorken-x, x 3> T/Q, or higher virtualities Q 3> Q s , 
the structure functions are exponentially small (since generated through tunelling). This 
strong suppression of the structure functions at large values of x and/or Q 2 implies the 
absence of point-like constituents in the strongly coupled plasma [ |32| , |33| , |26| , p7| . The critical 
value Qs{x) ~ T/x is known as the saturation momentum, since Eq. ( |3.21 ) is consistent with 



a parton picture in which partons occupy the phase space at Q < Q s (x) with occupation 
numbers of 0(1) p6[ . 

Returning to the flavor current of interest here, let us now identify the similarities and 
the differences with respect to the problem of the 7£-current, that we have just discussed. 

In the high-temperature case, where the tip of the D7-brane enters the BH horizon, there 
are no serious conceptual differences with respect to the 7£-current. For r>l, Eq. ( |3.15D is 
still valid, so the large-r dynamics is exactly the same as discussed in relation with Eq. ( |3.19D , 
At smaller r ~ 0(1), the EOM becomes more complicated (in particular because of the r— 
dependence of the profile function R v (r), which is non-trivial in that high-temperature case), 
but there is no ingredient in the dynamics which could prevent the fall of the flavor field 
into the BH. Hence, the appropriate boundary condition at p = 1 (the tip of the D7-brane) 
is still the infalling one, and moreover this condition can again be enforced ar large r> 1, 



where Eq. (3.19) applies. As before, this condition fixes c = i, thus finally yielding the same 



result for the DIS structure function as in Eq. ( 3.21 ), except for the overall normalization: 



F\ ~ ^ ^ f (—) (flavor current in the BH embedding) . (3.23) 



4r 2 (l/3) V6 

This is indeed the result found in 1 31 ] . In particular, the saturation momentum for the flavor 
current (in this high-temperature regime, at least) is exactly the same as for the 7£-current, 
cf. Eq. ( |3,22| ), since fully determined by the current interactions with the BH. 



7 One can show that Fl = i*2 — 1xF\ is parametrically of the same order as xF\ when x ~ x s , but it is 
relatively negligible when iCi s Efi]. 
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Consider now the low-temperature phase, which is the most interesting case for us here. 
For the DIS problem, it is natural to assume that Q > M gap S> T, or Q > Rq 1. The 
situation near the Minkowski boundary will be quite similar to that for the 7£-current: At 
relatively low energies Co <C Q 3 , there is a potential barrier at p > Q > Rq, which however 
disappears at larger energies Co > Q 3 - When this happens, the flavor field can penetrate all 
the way within the worldvolume of the D7-brane. However, this worldvolume ends up at 
p = Rq S> po, so there is clearly no possibility for this field to fall into the BH. Accordingly, 
the infalling boundary conditions do not apply here, but rather must be replaced by the 
condition of regularity at p = Rq (the tip of the D7-brane). In order to establish the fate of 



the high-energy current, we therefore cannot make the economy of actually solving Eq. ( 3.15 ) 
for all the values of p down to Rq. Still, there is an important 'technical' simplification which 
occurs at high energy: then, the potential barrier oc Q 2 plays no role anymore, so it is sufficient 



to consider the Q = version of Eq. ( |3. 15 ) . This equation also determines the spectrum of 



the light-like mesons in the plasma and, as we shall see, these two problems — the flavor DIS 
and the light-like mesons — are indeed strongly related to each other. 

4. Meson spectrum at low temperature 

In this section we shall construct piecewise approximations to the spectrum of the meson 
excitations in the low temperature phase, or Minkowski embedding, with the purpose of clar- 



ifying some global properties of the spectrum numerically obtained in Ref. [23] and exhibited 
in Fig. |l|. In particular, we shall follow the transition of the dispersion relation of a given 
mode from time-like to space-like with increasing momentum k, and thus identify the 'criti- 
cal' momentum k n at which the mode with radial quantum number n crosses the light-cone. 
Also we shall recover previous analytic results in the literature which concentrated on spe- 
cial limits, like zero-temperature [^3] or very high momentum [22|. The particular case of a 
light-like mode (co n (k) = k) will be given further attention in the next section, where we shall 
construct the exact respective solutions for both normalizable and non-normalizable modes, 
with the purposes of understanding DIS. 

4.1 Equation of motion in Schrodinger form 

For the subsequent analysis, it is convenient to change the definition of the radial coordinate 
once again, in such a way that the infinite interval < r < oo be mapped into the compact 
interval < ( < 1. Here £ is defined as 

^.2 ^,2 

£ = 9 1 \ r>2 i 9 ' (4-1) 

p z (r) RQ + r 2 

so in particular the Minkowski boundary corresponds to £ = 1. Then Eq. ( ^.15[) becomes 
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where we have replaced the name of the function by for more generality: indeed, Eqs. ( [3.15|) 
or apply not only to transverse vector mesons, but also to the pseudoscalar mesons 

corresponding to small fluctuations in the azimuthal angle <j) (the angle in the 89-plane 
transverse to the D7-brane; recall that the 'vacuum' embedding corresponds to eft = 0). 
Furthermore, we have defined 

where the virtuality Q 2 = k 2 — Co 2 can now take any sign (and thus the same is true for k 2 ). 
It is furthermore convenient to rewrite Eq. (|4.2| ) in the form of a Schrddinger equation, i.e., 
to remove the term involving the first derivative; this can be done by writing 

$(C) = iv(0- (4.4) 

The corresponding "Schrddinger equation" reads 

(-£>c + v £ (0)^(0 = 0, (4.5) 
1(1 + 2) , k 2 4n 2 (i-p 

m) = MRF + 4aw) r - ' (46) 

where we have allowed for one further generalization by adding to the potential the term 
corresponding to a generic value £, with I = 0, 1, 2, ... , for the angular 'quantum number' 
corresponding to rotations around the S^-sphere internal to the D7-brane. (Such rotations 
cannot be excited by either the flavor or the 7£-current considered in the previous section, 
so I = in the case of DIS. But modes with non-zero I can be excited by other operators 
in the boundary gauge theory, which are charged under the global SO(4) symmetry of the 



fundamental hypermultiplets |43|, 23, |25j.) Since the radial coordinate £ terminates at £ = 1, 
it is understood that the potential becomes an infinite wall at that point; given the structure 
of Eq. ( |4.6| ) , this additional constraint has no consequence except in the limiting case where 
t = and k 2 = 0. The zero temperature case is obtained by formally taking Q, = in 
Eq. (|J). 

Since we are interested in the normalizable modes describing mesons, we shall look for 
solutions ?pe(C) to Eq. Q4.5|) obeying the following boundary conditions [43] : 



1+2 

(— for C -> 

MC) oc { (4.7) 
(1 - C)~ for C 1. 

In what follows we would like to follow the change in the dispersion relation when increasing 
the meson momentum k for fixed quantum numbers (i.e., for a given mode). This study 
will drive us through different regimes in terms of the variables k 2 and Q 2 . The potential 
^(C) = V£=o(C) m these various regimes is illustrated in Figs. || and || 
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Figure 4: The potential V(() in Eq. (4.6) for £ = 0. Left: the case where fi <C k (as relevant at zero 
temperature, or for finite temperature but relatively small momenta k <C nR^). Right: the case where 
\k\ -C CI (this corresponds to modes with momenta k ~ ri-Rgj. 




Figure 5: The potential V(Q in Eq. ( |4.6| ) for £ = and relatively large, space—like, virtualities, with 
k < Left: k = 20. Right: k ~ 4£1 (tie limiting velocity regime). 

4.2 The low momentum regime: time like dispersion relation 

When the momentum k is sufficiently small (see Eq. ( 4.1 0| ) below for the precise condition), the 



dispersion relation is time-like (k 2 < 0) and it is such that the last term, proportional to Q 2 , 
in the potential becomes negligible (so that the potential has the symmetric shape shown in 
Fig. ^jleft). Then Eq. (42) is formally the same as at zero temperature and the corresponding 



solutions are known exactly [43]. Namely, the solution obeying the right boundary condition 



at £ = 0, cf. Eq. ( |4.7| ), reads (the overall normalization is chosen for convenience) 



t+-2 



MO = [C(l-C)]~ F(a,M + 2;C), (4.8) 
with a = i(2£ + 3+ \A + /U 2 ) , b = ^(2£ + 3- \A + /i 2 ) , 
where F(a, b; c; x) is the usual hyper geometric function, also denoted as 2-^1(0, b; c; x), which 



- 22 - 



obeys F(a,b;c;0) = 1 (see e.g. Chapter 15 in [65]) and we have set — k = /j, > 0. For 



Eq. (4.8) to also obey the correct boundary condition at £ = 1, the hypergeometric function 
must be regular at that point, which for the indicated values of the parameters a, b and c = £+2 
requires the hypergeometric series to terminate |5j|. That is, b = —n with n = 0, 1, 2, and 
then F(a, —n; c; x) is a polynomial of degree n in £. This condition yields the vacuum-like 
(i.e., T = 0) spectrum in Eq. ( 2.10Q , that is, 



it 



At finite temperature, Eq. (|4.9|) remains a good approximation so long as one can neglect the 
term oc il 2 in the potential, that is, for O 2 <C /^ 2 or, equivalently (recall (|4.3|) ), ui -C ni?j]. 
Since we also assume i?o ^ 1> it is clear that this condition is satisfied up to relatively high 
values of the momentum k, namely so long as 

k<nR%, or k <C n (M g 3 ap /T 2 ) . (4.10) 

Note that we include the radial quantum number within parametric estimates, e.g., fi n ~ n 
or ui n ~ ni?o, since we shall be also interested in large values n>l (whereas £ will never 
be too large). Although, for definiteness, we refer to the kinematical domain ( |4.10| ) as the 
'low-momentum regime', it is clear that, towards the upper end of this domain, the momenta 
are so large that k 3> M n ~ nM gap and thus the dispersion relation becomes nearly light-like. 

Still within this low-momentum regime, it is easy to see that the main effect of the 
"finite-T" term oc £l 2 in the potential ( |4.6| ) is to decrease the meson virtuality as compared 
to its "zero-T" value ( |4.9[ ). Indeed, a simple estimate for this effect is obtained by replacing 
fi 2 — > /i 2 + 16f2 2 in the l.h.s. of Eq. ( |4.9| ) (this procedure overestimates the correction when 
C — 1; but it should be correct at least qualitatively); hence the corrected virtuality reads 

{j?(k,n,£) w At( 0) (n,£) - 160^, with fj,f 0) (n,£) = 4(n + £+ l)(n + £ + 2) . (4.11) 

In particular, for non-relativistic momenta (k <C M n ~ nM gap ), the above dispersion relation 
can be expanded out as (the quantum numbers are kept implicit) 



k 2 A M (m 

u(k) « M rest + ^— , with M rcst « M (0) (l - 2/i^) , M kin « - _ ^ , (4.12) 



and therefore M res tMkj n ~ ^(o)- (-^(o) denotes the "zero-T" meson mass, as given by 
Eq. (2.10) or Q4.9| ).) The estimates ( 4.12j ) are in fact in agreement with the respective numer- 



ical findings in Ref. [g3| (see the discussion of Eq. (4.45) there). 

4.3 The intermediate momentum regime: light like dispersion relation 

With further increasing k, the energy u n of the mode n is also increasing and the last term, 
proportional to f2 2 , in the potential (|4.6| ) becomes more and more important. Since, at the 
same time, the virtuality // 2 of the mode is decreasing, it should be clear that for sufficiently 
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large k — namely, when k ~ uRq — one enters a regime where f2 2 S> // 2 and then the roles of 
the respective terms in the potential are interchanged: the term in Q 2 becomes the dominant 
one, while that in /i 2 represents only a small correction. Then the mode n is nearly light -like 
and in fact it crosses the light cone (i.e., its virtuality ^(k) changes sign) when varying k 
within this domain. We have not been able to analytically follow this transition, but the fact 
that it actually happens is quite obvious by inspection of the shape in the potential in this 
regime. This is shown in Fig. |4| right for the three cases of interest: (a) k 2 = —fj? is negative 
but small, (b) k 2 = 0, and (c) k 2 is positive but small. 

Namely, consider the genuine Schrodinger equation associated to this potential, that is, 

i-d 2 + v(o) v(o = mo , (4.i3) 

where E is the energy of a bound state. Given the shape of the potential in Fig. |] right, it 
is clear that bound states with both positive and negative energies will exist for all the three 
cases aforementioned. It is furthermore clear that, with increasing Q 2 at fixed k 2 the potential 
becomes more and more attractive, so some of the bound states will cross from positive to 
negative energies. This means that, for any fixed value of n 2 , there exist corresponding values 
of Q 2 such that the respective bound states have E = 0. These are, of course, the meson 
modes that we are interested in. 

This Schrodinger argument also suggest the use of the semi-classical WKB method for 
computing the meson spectrum. Given the shape of the potential this should be a reasonable 
approximation at least for sufficiently large numbers n>l (we set I = for simplicity). The 
Bohr-Sommerfeld quantization condition for the mode n with energy E n = reads 



nir 



Co 



d(V^ViC), (4.14) 



o 



where £o is the turning point in the potential, that is, (q = 1 when k 2 < and Co = 1 — 
when k 2 > 0. When k = 0, the integral is straightforward and yields 

O n w n, or ui n w V2uRq. (4.15) 

That is, the mode n crosses the light-cone at k = k n with k n ~ ^/2uRq. As we shall see in 
Sect, pi this is indeed the correct result when n 1. 

An interesting property of the spectrum near the light-cone, which will play an important 
role in our subsequent study of DIS (see Sect. |5|) and can be also understood on the basis of 
Eq. ( 4.14j ), is the extreme sensitivity of the dispersion relation to changes in the virtuality n 2 



around n = : so long as | k 2 \ <C fi 2 , a small change in k 2 entails a large change in ft 2 . Before 
we explain the origin of this property, let us first use Eq. ( [4.14 ) to render it more specific. 



Consider the space-like case k 2 > for definiteness, and denote e = «/4fl <C 1, so that the 
turning point lies at Cq = 1 — e. Changing the integration variable according to x = 1 — £, we 
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can successively write 



727T 

29 




a/x(1 - x) 



Q(x — e)^/x 2 ~ 



(4.16) 



7T 



,3/2 



dA 

71 



A - 0(A - A 2 - 1 



where we have observed that, after subtracting the dominant contribution tt/2 to the first 
integral, the subtracted integral is dominated by its lower limit x = e; this allowed us to 
perform the simplifications in the second line, where we denoted x = Xe. The final integral 
multiplying e 3 / 2 is clearly a positive number of 0(1). One can combine together the space-like 
and time-like cases into the following formula 



n 



9 i 



sgn(K 2 )CW — 



for 



< 9. 



(4.17) 



where C is a positive constant and sgn(x) = ®{x) — Q(—x) is the sign function. This formula 
shows that, when moving away from the light-cone, say, towards space-like virtualities, the 
energy of the mode grows by a substantial amount A9 n ~ 1 for a relatively modest increase 

1/3 

in the virtuality, from k = to k ~ 9n <C 9 n . In physical units, we change uj n by a large 
amount Au n ~ Rq 3> 1 when increasing Q from zero to Q ~ n 1//3 i?o ~ w^ 3 . 

This strong sensitivity of the dispersion relation to k? around k = can be traced back 
to the behavior of the potential near the turning point Co (recall that it is this turning point 



which controls the K-dependence of the integral in Eq. (4.14)). Namely, from Eq. ( |4.6j ) we 
deduce (for positive k with n <C O) 



dV 

8k 2 



C=Co 



n 

~ - > l 

K 



and 



d(-V) 



on 2 



C=Co 



(4.18) 



4.4 The high momentum regime: space like dispersion relation 



Consider now further increasing the momentum k of the mode n, beyond the critical value 
k n ~ \^2uRq at which the dispersion relation crosses the light -cone. Then the virtuality of the 
mode Q n will increase as well, i.e., the mode becomes more and more space-like (although, 
as we shall see, this virtuality remains relatively small, in the sense that Cu n ~ k ^> Q n ). For 
instance, Eq. (|4.17 ) implies that, so long as k <C 9, the virtuality grows with the energy (or 
the momentum) according to 



9 n -n 

Q"n. 



2/3 



or 



Qn 



p/3 
R 



\ 2/3 

1 



(4.19) 



Eventually, when 9 n 3> n, n n becomes comparable with Q. n and then the potential ( fOf ) has 
the shape shown in Fig. |5| (for t = and two different values of k). 
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As manifest on these pictures, when increasing the ratio k/Q, the turning point Co = 
1 — K./4Q in the potential moves towards £ = 0, i.e., towards the bottom of the D7-brane. 
Thus, clearly, the attractive region of the potential, which can support Schrodinger bound 
states with energy E n = (or, equivalently, space-like mesons), exists only so long as k < 40, 
and becomes very tiny (Co <C 1) when k approaches the upper limit 40 (cf. Fig. || right). What 
we would like to argue in what follows is that for sufficiently large momentum k S> uRq, the 
dispersion relation approaches this kinematical limit in which K n ~ 40 n : this is the 'limiting 
velocity' regime, previously mentioned in relation with Eq. ( 233 ) |23| , 22]. 



To that aim, let us compute the spectrum in the regime where k is indeed close to, but 
smaller than 40, in such a way that Co -C 1. The corresponding modes will be localized in 
the classically permitted region at ( < Co- It is then a good approximation to replace the 



d 2 c + £ -^-~^ + 2F 0. (1.20) 



potential (4.6) by its expansion near C = 0. The ensuing Schrodinger-like equation reads 

'1 + 2) _ 2e 
4C 2 C 

where we have denoted 

2e 2 = 4n 2 - — ~ m (n h 



4 V 4 

„2 



28 = 40 2 + — ~ 80 2 . (4.21) 

4 v ' 



When deriving Eq. ( |4.20| ) and also when simplifying the expressions in Eq. (4.21) we have 
anticipated the fact that, for the mode n, both $7 and k are very large but relatively close to 
each other, such that 

n n , Kn > n , n n - ^ ~ . (4.22) 



Also, we restricted ourselves to angular momenta £ O ri 

Eq. ( 4.20| ) is formally similar to the radial Schrodinger equation for a non-relativistic 



particle with mass m = 1 and electric charge e in the three-dimensional Coulomb potential 
Vc(C) = ( — e )/C; with — £ playing the role of the (negative) energy of a bound state. There 
are however some interesting differences with respect to the genuine Coulomb problem. First, 
the would-be 'angular' momentum of our fictitious 'Coulomb particle' is equal to 8 £/2, and 
hence it can also take half-integer values. Second, our radial variable £ is restricted to £ < 1 
and, moreover, the approximate equation ( |4.20| ) is valid only for £ <g 1; by contrast, in 
the corresponding Coulomb problem the radius £ can be arbitrarily large. Yet, this last 
difference should not be important for the situation at hand: given the potential barrier at 
C ^ Co) cf- Fig- !> it is clear that the actual meson wavefunction is exponentially decaying 
for £ > £o before exactly vanishing at £ = 1. When Co 1, there should be only a minor 



8 In the usual Coulomb problem in quantum mechanics, the centrifugal piece of the potential reads 1(1 + 1) /Q 2 



which is formally the same as that in Eq. (4.2C) provided one identifies I = 1/2 
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difference between the exact wavefunction, which is strictly zero at £ = 1, and its Coulombic 
approximation, which is exponentially small there. 



Hence, one can solve Eq. (4.20) by following the same steps as for the Coulomb problem 



in quantum mechanics [56]. The general solution which is regular at £ = reads 9 

l p e ( z ) = z ^ +1 e~ z/2 M(-v + l+ e/2,£ + 2;z), (4.23) 

z = 2V2£(, v = -i==. (4.24) 

v It 



Here M(a,b;z) is the confluent hypergeometric function, also denoted as \F\{a,b;z), which 
obeys M(a, b; 0) = 1 (see Chapter 13 in ||). Note that z ~ 4V§nC> cf. Eq. ( ^2l|) , hence C = 
1 corresponds to large z>l, where the asymptotic behaviour of Eq. (|4,23| ) becomes relevant. 
Specifically, for the solution to exponentially vanish at z > 1, the confluent hypergeometric 
series must terminate, which in turn requires 

-v + l + - = -n, n = 0,1,2,... (4.25) 

and then M(—n,£ + 2; z) is a polynomial in z of degree n (actually, a Laguerre polynomial 
up to a numerical factor M). The 'quantization' condition ( |4.25|) together with 



the definitions ( |4.21| ) and ( f4.24Q can now be used to deduce the meson spectrum in this high 



momentum regime. The resulting dispersion relation can be written in various, equivalent, 
ways, either function of the virtuality Q, 

Unl{Q) ~ ~ R 2 o Q ^ R% (n + 1/2 + 1) , (4.26) 
or as a function of the momentum k, 

2 ne (k) ~ vfk 2 * 8 (n + 1/2 + 1) A , (4.27) 

ito 

or, finally, in physical units: 

u? nl (k) - v 2 k 2 ~ (4^T) 2 (n + ^/2 + 1) -A_ . (4.28) 

-'"gap 

The last two equations feature the limiting velocity which has been generated here as 

^0 = 1-4, (4-29) 
K 



which is indeed consistent with Eq. ( [2.13 ) (recall that we assume Rq » 1). The results ([Dp 



( [4.28| ) are in agreement with a previous analytic study of this high-momentum regime, in 
Ref. p2|, which is more precise than ours. In any of these equations, the two terms appearing 



9 The notation z for the rescaled radial coordinate will be used only temporarily and should not be confused 
with the respective Minkowski coordinate, which is not anymore explicit in our analysis. 
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in the left hand side are large but comparable with each other, while the term in the right 
hand side, which expresses the deviation from the linear dispersion relation to = v$k and 
involves the dependence upon the quantum numbers, is comparatively small. 

By inspection of the meson wavefunction Eq. ( [4.23 ), where we recall that z ~ 4\/2^C) it 



is clear that the mode is localized near the bottom of the D7-brane, at C < l/^ n ^ 1- This 
domains lies within the classically allowed region at C < Co (indeed, Co = 1 — k/4£1 ~ n/tt n 
for a mode satisfying (|4.26| )), which confirms the consistency of our previous approximations. 
Interestingly, the higher the energy is, the stronger is the mode localized near C = (or 



p = Ro), in agreement with the numerical findings in Ref. [23]. 

Although the limiting velocity ( |4.29| ) is very close to 1 under the present assumptions, the 
virtuality Q 2 = k 2 — Q 2 of the meson is nevertheless very large, Q S> uRq, because its energy 
and momentum are even larger: uj n ~ k 3> nR\. Thus the mode looks nearly light-like in the 
sense that O n ~ k 3> Q, yet its virtuality is too high to be resonantly excited by an incoming 
space-like current: indeed, to be resonant with the meson, the flavor current should have an 
energy O and virtuality Q obeying u ~ (Rq/2)Q uRq, and therefore u>/Q 3 ~ (Ro/Q) 2 -C 1. 
According to the discussion in Sect. ||, such a current would encounter a large repulsive 
barrier near the Minkowski boundary and hence it would get stuck at large radial coordinates 
Pl^Q ^ n Ro, f ar away from the region at p ~ Rq where the would-be resonant mesons could 
exist. We conclude that such high-energy, space-like, meson excitations cannot contribute 
to the DIS of a flavor current. To investigate the possibility of DIS, we therefore turn to the 
only potentially favorable case, that of the 'nearly light-like' mesons with u n ~ k ~ uRq and 
arbitrarily small virtualities. 

5. Resonant deep inelastic scattering off the light— like mesons 

We now return to the problem of DIS off the strongly coupled J\f = 2 plasma at low temper- 
ature, as formulated in Sect. ||. Recall that we are interested in a relatively hard space-like 
flavor current, with virtuality Q > Rq (or k > 1). So long as the energy ui of this current 
is relatively low, ui <S Q 3 (or f2 <C k 3 ), there is a repulsive barrier which confines the dual 
gauge field Ai(p) near the Minkowski boundary, where no resonant meson states can exist. 
(This barrier is visible in Fig. ||| as the repulsive potential at ( > Co-) We thus conclude that 
the flavor structure functions vanish when lj <^ Q s , so like for the 7£-current. 

However, the situation changes when the energy of the current is sufficiently high, such 
that Co > Q 3 (or Q > k 3 ). Then, the repulsive barrier becomes so narrow that it plays no role 
anymore (this is visible as the curve 'k 2 > 0' in Fig. |] right). Indeed, even in the presence of 
this barrier, the gauge field can penetrate across the barrier, via tunneling, up to a distance 
orl — C~1/k 2 ; when Q > k 3 , this penetration is larger then the width 1 — Co = k/AQ 
of the potential barrier, and then the field can escape in the classical allowed region at C < Co- 
Thus, the field has now the capability to excite vector mesons at any value of p within the 
wordvolume of the D7-brane. Moreover the kinematics of the high-energy current matches 



with that of the nearly light-like mesons discussed in Sect. 4.3. Indeed, those mesons have 
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energies u n ~ nR,Q, cf. Eq. ( 4.15 ), which can match the energy (D > (J 3 > of the current 
with a suitable choice for n. Furthermore, for a given n, there are mesons at all virtualities 
Qn ^ nRoi and in particular such that Q 3 < u) n , so like for the current. 

These kinematical arguments indicate that the high-energy flavor current can disappear 
into the plasma by resonantly exciting nearly light-like mesons. In what follows, we shall 
demonstrate that this picture is indeed correct, by explicitly computing the decay rate (i.e., 
the imaginary part of the current-current correlator) corresponding to the resonant excitation 
of large-n light-like mesons. The restriction to large quantum numbers n> 1, that is, to very 
high energies u> 3> Rq, is necessary for technical convenience, but it also has the advantage 
to make the physics sharper. In that case, it becomes possible to smear our the delta-like 
resonances associated with the individual mesons and thus obtain a spectral function which is 
a continuous function of a) and hence describes DIS. Remarkably, that spectral function turns 
out to be identical with the DIS structure function in the high-temperature phase ('black hole 
embedding'), Eq. ( ggp , that was previously computed |3l|] by imposing infalling boundary 
conditions at large p 3> i?o (cf. the discussion in Sect. [T^). 

5.1 Light like mesons: exact solutions 

In this subsection we shall concentrate on the EOM for light-like, transverse, gauge fields 
in the worldvolume of the D7-brane, that is, Eq. ( |3.15 ) with Q = 0, for which we shall 



construct exact solutions obeying the condition of regularity at r = (or p = Rq). By using 
the asymptotic expansion of these solutions at large p and high energy, we shall study their 
behaviour near the Minkwoski boundary and thus distinguish between normalizable and non— 
normalizable modes. In particular, this procedure will yield the spectrum of the light-like 
mesons for large quantum numbers n> 1. 



Once again, it is more convenient to use the variable £ defined in Eq. (4.1) and which 



has a compact support. Then the relevant EOM is Eq. (O) with k = 0, or, equivalently, 
the 'Schrodinger equation' (|4.5| ) with £ = and k = 0. We shall choose the latter, that we 
rewrite here for convenience: 



di- 4 ^ C) )V(C) = 0- (5.1) 



)2 



Clearly, this is a particular case 10 of Eq. ( |4.20| ) that we have solved already, namely it is the 
limit of that equation when I = and 2e 2 = 2£ = AO, 2 . The solution which is regular at 



C = is then obtained by adapting Eq. ( 4.25 ), and reads 

<J?(C) = J V(C) = e" 2C? M(-n + 1, 2; 4^) , (5.2) 

which once again features the confluent hyper geometric function M(a,b;z). This solution 
takes on a finite value on the Minkowski boundary at £ = 1 and in general, i.e., for generic 



In spite of this formal similarity, one should keep in mind that Eq. (4.20) and, respectively, Eq. ( |5.l[ ) apply 



to different physical regimes. In particular, Eq. (5.1) holds for any £ < 1, whereas Eq. (4.2C) is valid only for 

C«i. 
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values of the energy parameter 17, it represents a non-normalizable mode. The normalizable 
modes describing mesons are obtained by requiring the solution to vanish at £ = 1: 



M(-17 n + l 5 2;417 n ) = 0, (5.3) 

with 17 n the on-shell energy of a light-like meson. For generic values 17 ~ 0(1), this equation 
is difficult to solve except through numerical methods. A similar mathematical difficulty arises 
when trying to use Eq. ( |5.2[) in order to compute the current-current correlator according to 
Eq. ( grp . For all such purposes, one needs the behavior of the solution near the Minkowski 
boundary at £ = 1, and this is generally difficult to extract from Eq. ( f>.2\) . 

However this mathematical problem becomes tractable for the high energy regime of 
interest here, which is such that Q 3> Rq, or 17 3> 1. Then, one can use a special asymptotic 
expansion of the function M(a,b;z) with a < 0, which applies when the variables \a\ and z 
are simultaneously large and such that z ~ 26 — 4a 3> 1 . This last condition is truly essential, 
since in general, i.e., for generic values of \a\ and z which are both large but uncorrelated with 
each other, very little is known about the asymptotic behavior of M(a,b;z). This specific 
limit is precisely the one that we need for our present purposes: indeed, in Eq. (|5.2| ), we have 
z = 417£ an d 26 — 4a = 417, and therefore z ~ 26 — 4a = 417 3> 1 in the high energy limit and 
in the vicinity of Q = 1. The asymptotic formula which applies to this case is formula 13.5.19 
in Ref. l55|| and can be formulated as follows: when 



(26 - 4a) 



i + J 



with |f | < (6 -2a) 2 / 3 , (5.4) 



(6 -2a) 2 / 3 . 

then (below, Tip) is the Euler function, and Ai and Bi are the Airy functions) 

M(a,b;z) = el (6 - 2a) ^ b T(6) [Ai(£) cos(vra) + Bi(£) sin(vra) + 0(|6-2ar 1 )] . (5.5) 
In order to adapt this formula to Eq. (|5.2|), we shall write 



so that the variable f be positive. Then for 17 3> 1 and 1 — £ <C 1, the solution ( p.2[ ) becomes 
(up to an irrelevant overall normalization) 

$(0 ~ Ai(-f) cos(«tt) - Bi(-f ) sin(nTr) for < f < (217) 2 / 3 . (5.7) 

In particular, for relatively large ^ > 1, one can use the asymptotic expansions of the Airy 
functions to deduce 

<&(£) ~ -^-^ sin {^e 12 + ~ - «tt) for 1 « f « (217) 2 / 3 . (5.8) 



The solution ( p.7|) has the same general structure and validity range as the approximate 
solution shown in Eq. ( |3.19| ) — in particular, the argument f of the Airy functions is indeed 
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the same in both equations, as it can be checked by using Eqs. (JO]) and (|5.6D — , and this 
should not be a surprise: as explained in Sect. ||, Eq. ( |3.19| ) is the general form of the solution 
at high energy and large p. The whole purpose of a more complete analysis at smaller values 
of p, like the one that we have just performed here, is to fix the coefficients of the two 
Airy functions appearing in that equation. Note that, unlike for the solution with infalling 
boundary condition, i.e., Eq. ( p. 19 ) with c = i, the coefficients in Eq. ( |5.7| ) depend upon the 



energy variable O. 

As a first application of Eq. (|5.7|) , we now use it to determine the energies of the light -like 
meson excitations according to Eq. (fO|). To that aim, we also need fl55fl 



Ai(0) = „,„ 1 ; : = -= Bi(0) , Ai'(O) = — T7 — — = — 7= Bi'(O) , (5.9) 
v ; 3 2 / 3 r(2/3) V3 3V3r(l/3) ^3 

(the formulae involving the derivatives will be useful later on). Then a simple calculation 
shows that for £ = the right hand side of Eq. ( |57^ ) is proportional to sin[(f2 — l/6)7r], and 
then Eq. ( ^^ ) immediately implies 

n n = n+l, (5.10) 

O 



in agreement with Eq. ( 4.15 ). It can be numerically checked that Eq. ( 5.10 ) is a very good 
approximation to the zeroes of Eq. ( |5,3D already for small values n = 1, 2, .... One may think 
that the constant shift O n — n = 1/6 in the eigenvalues is merely a tiny correction that can 
be safely ignored at large n, but this is generally not the case. Note first that this shift is 
uniquely determined by the values of the two Airy functions at £ = 0, as it can be checked by 
inspection of the previous manipulations, and hence it is not affected by the approximation in 
Eq. ( |5.7| ) . Moreover it is essential to take this shift into account whenever one is interested in 
the behavior of the solution near O = D, n , which will be also our case in the next subsection. 
The wavefunction corresponding to the mode n is obtained by replacing — ► O n within 



the general formulas (|5.2|) or fl5.7|). For radial coordinates deeply inside the D7-brane, where 



the asymptotic expansion fl5.7| ) does not apply, one can rely on the exact solution (5.2), but 
this is perhaps a little opaque. An approximate formula valid for intermediate values of £ will 
be constructed via the WKB method in Appendix [A|. This WKB solution, shown in Eq. (A.2), 



is consistent with the asymptotic behaviour (5.8) and has the nice feature to exhibit exactly 



n nodes in the interval < ( < 1, as a priori expected for the nth radial excitation. 
5.2 Current current correlator and DIS 

We are now prepared to compute the current -current correlator for a highly energetic, nearly 
light-like, current, and thus make the connection to DIS, as anticipated. To that aim, we shall 
use Eq. ( 3.12j ) together with the asymptotic expansion of the solution near the Minkowski 



boundary, Eq. (|5.7| ). Specifically, Eq. (|3.12[) involves (recall that Ai = in our present 



notations, and the radial coordinate r in Eq. ( 3.12 ) is related to the variable £ which appears 
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in Eq. flO|) via Eqs. and (|5j 



] -- >(ir 






J 1 — >oo 







^3 d r Aj(r,uj,k ) 
Ai(r,uj,k) 

Then a straightforward calculation using <!>(£, 0) from 

3nN c N { T 2 



ni(w) 



2r 2 (l/3) v6/ 



— = + COt ( I 7T 



(5.11) 

together with Eq. yields 

(5.12) 



with Q related to Q via Eq. (|4.3| ). As expected, the function Ili(d)) exhibits poles at the 
energies O n = n + 1/6 of the light -like meson modes. These poles can be made more explicit 
by using the expansion of the cotangent as a series of simple functions: 



cot TTX 



1 2x^ 1 
— + — > 

TTX 7T L — ' 3T — 

n=l 



(5.13) 



To extract the spectral weight associated with these poles, i.e., the imaginary part of the 
correlator, we use retarded boundary conditions, to — > uj + ie, together with the formula 



lim 



1 



n 



V — — iir5(x 
x 



n) 



+o x + ie 

We thus find (recall that our energy variable is always positive) 

2/3 00 



37T 

u 7 2r 2 (i/3) 



(I)' 



^ 5(0 - 0„) . 
=i 



(5.14) 



(5.15) 



This result has been obtained here by working with a light-like current, but a similar result 
holds also for a highly energetic space-like current with u> > Q 3 , since the repulsive barrier 
plays no role in that case, and since the plasma can indeed sustain slightly space-like mesons 
which are resonant with the current. (In fact, the WKB method in Appendix |A] can be easily 
generalized to such slightly space-like mesons.) The emergence of the delta-functions in the 
imaginary part of the current -current correlator confirms our expectation that a high-energy 
flavor current can resonantly produce mesons in highly excited states (n S> 1) and thus 
disappear into the plasma. 

Taken literally, Eq. ( [5.15| ) would imply that the meson production by the current can 
only occur for a discrete set of energies which are resonant with the energies uj n = 
l)Rl of the light-like meson excitations in the plasma. However, an argument based on the 
uncertainty principle together with the peculiar structure of the meson dispersion relation 
near the light-cone (cf. Sect. |4.3j ) shows that in order for the absorbtion process to be 
observable, one needs to average Eq. ( |5.15| ) over neighboring levels. 

The argument goes as follows: A current with a given, large, momentum k which is 
produced by a source acting over a finite time interval 5t has an uncertainty 5uj ~ 1/St 
in its energy, and hence an uncertainty 5Q 2 ~ 2k5uj in its virtuality. (We have used here 
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k = \/ u> 2 + Q 2 ~ uj + Q 2 /2u) at high energy k S> Q.) As we shall demonstrate in Sect. ||, via 
an analysis of the time scales for the current interactions in the plasma, the typical interaction 
time for a nearly light -like current scales with its momentum like 

tint ~ ^ V3 • (5.16) 

(As usual, a bar over a kinematic variable denotes the dimensionless version of that variable 
measured in units of ttT, e.g. i = nTt.) So, for this process to be experimentally observable, 
the source producing the current must act over a comparatively short period of time: St < ti nt . 
(If 5t 3> ti n t , one cannot distinguish between the absorbtion of the photons in the plasma and 
their reabsorbtion by the source.) This in turn implies 

Su ~ J- ~ £173 ' or ^ > fc 1/3 • (5-17) 

This means that the high-energy flavor current has the potential to produce meson excitations 
with momentum k (the momentum of the current) and virtualities within a range 5Q > fc 1 / 3 
around Q = 0. 



At this point one should remember the discussion towards the end of Sect. 4.3, about the 
high sensitivity of the nearly light -like meson dispersion relation to changes in the virtuality. 
Let us rephrase that discussion but from a different perspective: assume that the momentum 
k of the meson is now fixed, but consider changes in the mode quantum number n associated 
with changes in virtuality k (near k = 0). As it should be clear from Eq. (4.17), n varies by 



a number of order one when k changes by Ak ~ fi 1 / 3 , that is, when Q varies by AQ ~ /c 1 / 3 . 
This is of the same order as the lower limit on the uncertainty ( |5.17|) in the virtuality of the 
current. Thus, for a given momentum k, the current has the possibility to be resonant with 
several meson states, with neighboring quantum numbers. 

Thus, in order to compute the total interaction rate for the current, as given by the 
imaginary part of the current-current correlator, we are allowed to average Eq. ( gig ) over 



several neighboring levels and thus smear our the delta-function resonances. This averaging 
amounts to integrating Eq. ( |5.15 ) over an interval 5Q which contains a few neighboring 



resonances and dividing the result by 5Q, thus yielding the following result 

ImI,i « = ir*m NMT U) • <5 - 18) 

for the imaginary part of the retarded 2-point function of a high energy flavor current with 
momentum k and energy u) ~ k. In particular, when the current is space-like (with relatively 
small virtuality, though: Q < A; 1 / 3 ), a non-zero imaginary part is synonymous of deep inelastic 
scattering, and the above result can be identified with the DIS structure function: F\ = 
(l/27r)lmlli. By comparing this result to Eq. ( |5.15| ) we see that the structure function 
in this low temperature phase is exactly the same as in the high-temperature phase, or 
'black hole embedding'. This coincidence reflects that fact that in both cases the current is 
completely absorbed into the plasma, although the respective mechanisms are quite different: 



- 33 - 



resonant excitation of nearly light -like mesons at low temperature and, respectively, partonic 
fluctuations (a quark-antiquark pair together with arbitrary many M = 4 quanta), which 
disappear into the plasma via successive branching, at high temperature. This physical 
picture will be further clarified by a discussion of the relevant time scales in the next section. 

Although our present calculations apply to the transverse field and structure function 
alone (recall that Fx = 2xFi), it is clear that a similar argument must be valid also in the 
longitudinal sector. So the corresponding, flavor, structure function Fl = F2 — Ft can be 



deduced by simply rescaling, by a factor 4Nf/N c , the respective result for the 7\!.-current 26] 



N c N f Q 2 x 

Fl = F2 — 2xF, = . (5.19) 

To summarize, for the N = 2 plasma at strong coupling, the above results for the flavor DIS 
structure functions are valid at either low, or high, temperatures, for high enough momenta 
k S> Rq and for sufficiently low virtualities Q < k 1 ^ 3 . In physical units, these conditions 
amount to k 3> T(M gap /T) 3 and Q < Q s (x), where the saturation momentum Q s {x) ~ T/x 
is the same as for the 7£-current. 

6. Time dependence and physical picture 

In what follows we would like to provide a space-time picture for the interactions of the 
flavor current in the strongly-coupled N = 2 plasma, and thus in particular clarify the 
energy averaging over neighboring resonances performed in the previous section. To that 
aim, we need to assume that the source producing the current has acted over a finite interval 
of time 8t, which is much shorter than the typical interaction time in the plasma, i; nt , that 
we shall compute. Accordingly, the current is not a simple plane-wave anymore, but rather 
a wave-packet in energy. To study the dynamics of this wave-packet, it is convenient to first 
reformulate the respective EOM as a time-dependent Schrodinger equation. This will also 
give us insight into the typical time scales for meson excitations in the plasma. 

A current produced over a finite period of time 5t can be described as a wave packet in 
energy, with a width 5oj ~ 1/St which is much smaller than the central value loq ~ k. (We 
assume that the current has a sharply defined longitudinal momentum k and we consider the 
high energy kinematics where to ~ k.) To study the evolution of this wave packet with time, 
we need to restore the time-dependence in the respective equations of motion, as written 



down in Sect. 3.1. For the problem at hand, this can be readily done by replacing 

u? -> (oj + id t ) 2 ~ loq + 2iLO d t , (6.1) 

in equations like (|3.15| ). Indeed, for the (relatively narrow) wave packet under consideration, 
the time dependence will be represented by a wave e~ luJ(>t modulated by a relatively slowly 
varying function (\idt\ <C ujq). Starting with ( |3.15 ), we thus obtain 



- —r dfAi = Ai + - At + — \ + )Ai , (6.2) 
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where we recall that t = irTt, k = k/irT, p 2 = Rq + r 2 , and a dot denotes a derivative w.r.t. 
r. In writing this equation, we have replaced ujq by k everywhere except in the virtuality 
term Q 2 = k — ujQ (which in what follows will be treated as a small correction) . Also, we have 
neglected a subleading term involving the time derivative which is proportional to 1/p 8 (recall 
that p > Rq ^> 1). We would like to rewrite this equation as a time-dependent Schrodinger 
equation, since then we can rely on the techniques and intuition developed within quantum 
mechanics. The first step is to eliminate the term oc A{, by writing Ai(i,r) = &(t, r)/r 3 / 2 : 

^lfi~\W 2 + ^P 2+ ^ r ~^F) ( ' 

Given the r-dependent factor l/p 4 multiplying the time derivative in the l.h.s., this equation 
is not yet in Schrodinger form. It turns out that the canonical, Schrodinger, form of the 
equation can be achieved by changing the radial coordinate one more time, namely, by using 
the angle 6 introduced in Sect. 0, cf. Eq. (|2.7| ), to that purpose. Specifically, after writing 



t = 2V2R t , p = , $(t, r) = ^— (i, 8) , (6.4) 

sm sm 

one finds that Eq. ( |6.3| ) takes the canonical form (with £1 and k defined as in Eq. ( |4.3|) ) 

f I d 2 \ 

n#) = _^_ 80 ^, + ^_ ? jL_ + ^. (6 . 6 ) 

Formally, these equations describe the quantum dynamics of a non-relativistic particle with 
mass O and Hamiltonian H defined by the r.h.s. of Eq. ( |6.5| ). 

Note that the term involving the virtuality k within the potential is independent of 8 
and thus merely acts as a constant shift in the total energy, in the same way as the time 
derivative. This is as expected: the time-dependence in the problem arises because of the 
uncertainty 5uj in the energy of the wave-packet; for a fixed momentum k, this corresponds 
to an uncertainty 5Q in the virtuality, such that 5oj ~ 5Q 2 /2k. In the units of Eq. ( |6 . 5|) , this 
amounts to id T ~ 5k 2 /2fi, so the time-derivative and the (central value of the) virtuality 
act indeed on the same footing. This being said, in what follows we shall often ignore the 
last term k 2 /2Q in Eq. ( |6.6[) , precisely because we are interested in situations where the 
fluctuations in virtuality associated with the uncertainty principle are larger than the central 
value Q 2 . That is, we shall assume Q <€. k 1 / 3 (the high -energy kinematics where the DIS 
process becomes possible), whereas we shall see that 5Q > k 1 ^ 3 . 

The potential in Eq. ( |6.6| ) is displayed in Fig. |6] for the interesting case f2 3> 1 and tz = 0. 
The Minkowski boundary lies at 9 = and the bottom of the D7-brane at 6 = ir/2. Other 
remarkable points are the two classical turning points, 0\ and 63, and the minimum of the 
potential at 82- One finds 

3V6 vr (3/2) 1 /4 31/2 

^^W, ^-0 3 ^—, (6.7) 



2nV3 ' 2 (sny/ 2 2 sr> 



- 35 - 




Figure 6: The potential V{9) in Eq. (3.6) for Q — 5 and k — 



so that 9\ <C 1, while 9 2 and #3 are close to tt/2. The potential has a rather deep minimum: 
V{9 2 ) ~ -80. 

To estimate the interaction time for the flavor current, we shall evaluate the time that 
the wave packet takes to travel from the boundary to the interior of the D7-brane, where it 
can excite mesons. A similar analysis for the case of the 7^-current was presented in Ref. |2j 
Near the boundary, where 9 < 1, the equation becomes 

.d^ ( 1 d 2 3 1 , , 

(-^t^ + t^K (6-8) 



dr V 20 d9 2 80 9 2 
with the following, exact, solution (C is a constant) 

¥(r,0) = C — expfi— J , (6.9) 

which describes diffusion: at early times, the penetration 9 of the wave packet in the radial 
dimension grows with r like 

9(r) ~ for r < n ~ O 1 / 3 . (6.10) 



Here t\ is the time the wave takes to reach the point 0\, cf. Eq. ( |6.7D , where the potential 
becomes attractive; in physical units, t\ ~ i?oTi ~ A; 1 / 3 . 

For t > T\, and so long as 9 lies in between 9\ and the minimum of the potential at 9 2 , 
the wave packet falls in the potential, essentially by following the classical equation of motion: 
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For most of this travel, 6 is still small, so sin 4 ~ # 4 and 

d 2 9 „, , 1 



^" 32g J =>^)" i/ gl - 4 (r-n) - (6 - 12) 
So, the penetration becomes of 0(1) at t ~ T2, with 

r 2 - Tl ~ -i- ~ ft 1 / 3 => t 2 - Fl ~ h ~ fc 1 ^ . (6.13) 

4t71 

That is, it takes (parametrically) as much time to the wave to fall in the potential up to 
large distances 6 ~ 0(1) as it takes to diffusively reach the attractive part of the potential, 
although the respective distances are very different: 82 3> 8\. This is so because along the 
second part of the trajectory, from 8\ to 82, the wave has un accelerated motion under the 
influence of the BH. 

This time i 2 ~ k 1 ^ is already a realistic estimate for the interaction time, since the 
current can resonantly produce mesons when 8 is of 0(1). Moreover, this estimate would not 
change if the meson was to be produced further down in the worldvolume of the D7-brane 
(say, near its lower tip at 8 ~ tt/2), since the final part of the fall in the potential is very 
rapid because the potential is so attractive around 62- We thus conclude that t- mt ~ fc 1//3 , 



as anticipated in Eq. (5.16|). This also confirms that the typical uncertainties in the energy 



and the virtuality of the wave packet satisfy 5ui > 1/fc 1 / 3 and, respectively, 5Q > fc 1 / 3 , in 
agreement with our original assumptions. 

Since determined by the dynamics relatively close to the Minkowski boundary, this inter- 
action time is independent of Rq, and thus is parametrically the same as for the 7£-current 
p7f |. As we shall shortly see, a similar conclusion holds also for the light-like meson excita- 
tions with large quantum numbers n ^> 1 : such a meson has a period (the interval of time 

- 1/3 

corresponding to one rotation around a semiclassical orbit) Ai n ~ Un and spends most of 
this time at radial locations far away from Rq, namely around p ~ u^^Rq. The subsequent 
calculation will also shed light on an interesting property of the meson spectrum, which played 
an important role in our previous argument: the strong sensitivity of the dispersion relation 
to variations in the virtuality Q around the light-cone (Q = 0). 

To that aim, we resort again on the Bohr-Sommerfeld quantization formula, valid for 
large n. We thus write 

nir = / d6 P {9-E n ,Q), p(8;E,n) = y/2Q(E - V(Q)) , (6.14) 



where E n are the energy levels associated with the Schrodinger Hamiltonian in Eq. (6.5). 
These energies are defined by the usual eigenvalue problem H^ n = E n *$> n and depend upon 
the two parameters and k within H. They should not be confused with the meson en- 
ergies Q n , which rather correspond to the special values of f2 (at a given k) for which the 
homogeneous equation = has non-trivial, normalizable, solutions. One clearly has 

E n (n,K)=0 for n = n n (n). (6.15) 



-37- 



This last equation has a unique solution for any n = 0, 1, 2, as it can be easily recognized 
by inspection of the potential in Eq. ( |6.6[ ) and Fig. ^. By choosing E n = in Eq. ( |6.14 ), 



one could work out the corresponding integral and thus recover our previous result that, e.g., 
^n(^) = n, cf. Eq. ( 4.15| ). However, for the present purposes, we shall rather use Eq. ( |6.14 ) 



to compute the level spacing of the Schrodinger energies near E = 0, that is 

AE(n n ) = E n+1 (n n ) - E n (n n ) = E n (n n ) ( K = o), (e.ie) 



Using Eq. (6.14), this is obtained as (we anticipate that AE <C 1) 

vr = / de[p(e;AE,n n )-p(6;0,n n )] ~ AE / d9 -§ = AE d9 — (6.17) 

Indeed, this quantity AE(tt n ) provides the answer to the two questions that we are interested 
in, as we explain now: 

(i) The last integral in Eq. (|6.17| ) is the same as half of the period for a round trip 
around a semiclassical orbit: indeed, p = £ld T 9, hence (Q/p)d9 = dr. 

(ii) The quantity AE(£l n ) characterizes the response of the meson dispersion relation 
to variations in the meson virtuality near the light -cone. This can be understood by recalling 
the discussion below Eq. ( |6.6D , about the last term, k 2 /2J7, in the potential. More precisely, 
AE(Q n ) is a measure of the increase in the virtuality which is needed to jump from one meson 
level (n + 1) to the neighboring one (n) at fixed energy and in the vicinity of the light -cone. 
In formulas, AE(Q n ) = K^/2Q, n with K n defined by Q n+ i(0) = Q n (K n ), or 

E n (n n (K n ),K n ) = E n+1 (Q n (K n ),0) = 0. (6.18) 

Returning to Eq. Q6.17| ), this can be evaluated as 

- i r0 3 / Q i \ -I/ 2 

JL ~ i / d6 sin 4 6 - — -r 9 (6.19) 
AE 4./fl, V MQ^ sin 2 9 cos 2 9 J V ; 



n 

d99 1 



4j 9l ^W^Of 24^ 



B(l/6,l/2), 



where in the second line we have used the fact that the integral is dominated by its lower 
limit 9± <^ 1, and B(l/6, 1/2) is the respective Beta function. 

We thus find that the period for motion of a meson around the semiclassical orbit with 

_ y /3 

quantum number n is At n ~ \/9\ ~ u> n ■ Moreover, the fact that the above integral is 
dominated by 9 ~ 9\ also means that a quantum particle in the bound state with energy 
E n ~ spends most of its time at relatively large radial distances p ~ u^ 3 ~ n l / 3 Ro- This 
implies a similar property for the light -like meson with energy Co n ~ \^2nR.Q. Of course, the 
meson wavefunction has support everywhere in the range Rq < p < ii^^Rq, but the radial 
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velocity d T 9 is smaller towards the upper end of this range (as it should be clear by inspection 
of the potential in Fig. |6|), therefore there is a larger probability to find the meson in that 
region than towards the bottom of the D7-brane. 



Furthermore, the result for AE in Eq. (|6.19[) implies 



2Q r 



n 2 J 3 



(6.20) 



1/3 



This is in agreement with our previous discussion of Eq. (4.17) in Sect. \L 



or Q Ti 

there is enough to make a rather small change AQ n ~ n 1/3 Ro in the meson virtuality in order 
to jump from one mode to another at fixed energy, whereas one needs a substantially larger 
increase in the energy of the meson, namely ACu n ~ B^, in order to make that jump at fixed 
virtuality. The above result for Q n is moreover of the same order as the fluctuations 5Q in 
the virtuality of the flavor current due to its energy uncertainty, thus justifying the energy 
averaging performed in our previous calculation of the spectral weight, in Eq. ( 5.1§p . 
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A. WKB approach to light— like mesons 



In this appendix, we shall use the WKB approximation to rederive the results in Sect. 5.1 on 



the spectrum and the wavefunctions of the light -like mesons in the limit of a large quantum 



number n ^> 1. We have already noticed in Sect. 4.3 that the WKB method correctly 
reproduces the spectrum fl5.10| ) except for the constant shift 1/6 in the eigenvalues Q n (which 
is of course a subdominant term at large n). As we shall see, this constant shift is also the 
only source of deviation between the WKB eigenfunctions and the asymptotic expansion of 
the corresponding exact result, as given by Eq. (| 



Specifically, let us view Eq. (5.1) as the Schrodinger equation for a non-relativistic quan- 
tum mechanical particle in a stationary state with energy E = 0. We recall that the potential 
in that equation becomes infinite at the end point £ = 1 (cf. Fig. |] right). Then for Q < 1, 
the associated wavefunction can be obtained in the WKB approximation as [^] 

V>(C) = ~J= sm (jf dC'KO + a) , p(C) = v^V(Cj . (A.l) 
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The shift a in the argument of the sine function is not accurately determined by the WKB 
method, but will be later fixed by matching onto the exact solution near Q = 1. With the 
potential V(() in Eq. (|5.1[) , the integral in Eq. ( A.l ) is straightforward and yields 

1/4 

Vc- 



m = c 



c 



i-C 



sin 



20 



arcsm 



Vc(i-c) + a 



(A.2) 



This solution is not reliable very close to the end points at Q = (where the potential becomes 
singular) and Q = 1 (where the potential is not differentiable) , but it should be a reasonable 
approximation at the intermediate points. A similar WKB solution can be constructed on the 
right of C = 0. The condition that the two solutions match with each other at intermediate 
points leads to the Bohr-Sommerfeld quantization condition ( [4.14 ), which in turn implies 
O n ~ n when ti>l. With O n = n, the WKB solution (A.2) has exactly n nodes in between 
and 1. 

When approaching the end point at ( = 1, we expect Eq. ( A.2| ) to remain a good ap- 
proximation so long as 56] 

1 __ 1 
" 20' 



(1 - C) 3/2 » 



(A.3) 



where —dV/dC, = 40 2 is the left derivative of the potential at £ = 1. This range of validity, 
which in terms of the variable £ introduced in Eq. fl5.6j ) amounts to ^ 3> 1, is wide enough 
to allow for a matching between the WKB solution and the exact solution near £ = 1. The 
latter is the following linear combination of Airy functions 



MO * ^ Ai (-0 - ~Bi( 



-£) for < £ < (20) 



2/3 



(A.4) 



where the relative coefficient has been fixed by the condition that ■0 = when £ = 0. Clearly, 
Eqs. (A.2) and (A.4) have a common validity domain at 1 -C £ ^ (2Q) 2 / 3 , and within that 
domain they are indeed consistent with each other, as we now check. 

To that aim, we need to compare the approximate form of Eq. ( |A.4| ) valid at £ 3> 1, 
which is obtained similarly to Eq. (|5.8[) and reads 



V>b(£) ^ 



C 
^74 



sin 



(A.5) 



to the approximate form of Eq. (A.2) valid near £ 
obtained as 



1 (i.e., for £ <ti (20) 2 / 3 ), which is readily 



C 
^T74 



sm 



e /2 + a) . 



(A.6) 



Clearly, these two expressions, ( |A.5| ) and (A.6), are consistent with each other, as anticipated. 
Moreover they perfectly match provided one choses a = vr/12, thus completely fixing the 
WKB solution ( A.2 ), valid for any £ which obeys the condition ( |A.3| ) and which is not too 
close to C = 0. 
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